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Menetelma vastaanottimen tahdistamiseksi ja vastaanotin 

Nyt esilla oleva keksinto kohdistuu oheisen patenttivaatimuksen 1 joh- 
danto-osan mukaiseen menetelmaan vastaanottimen tahdistamiseksi, 
5 oheisen patenttivaatimuksen 17 johdanto-osan mukaiseen vas- 
taanottimeen seka oheisen patenttivaatimuksen 32 johdanto-osan 
mukaiseen elektroniikkalaitteeseen. 

Eras tunnettu paikanriusjarjestelrtia on GPS-jarjestelma (Global 
10 Positioning System), joka kasittaa yli 30 satellllttia, joista samanaikal- 
sesti vastaanottimen nakyviss^ on maksimissaan 12. Nama sateliiitit 
lahettavat mm. satelllitin ratatletoa (Ephemeris data) seka tletoa 
satelliitin kellonajasta, Sijainnin maarityksessa kaytettava vastaanotin 
.pad.ttelee sijaintinsa normaalisti siten, etta vastaanottimessa lasketaan 
15 useammasta paikannusjarjestelmaan kuuluvasta satelliltista samanai- 
kaisesti lahetett^van signaalin kulkualka vastaanottlmeen. Sijainnin 
mddrittamiseksi on vastaanottimen vastaanotettava tyypillisesti vahin- 
taan neljan n^kyvissa olevan satelliitin signaali, jotta stjainti voidaan 
laskea. 

20 

Jokainen GPS-jarjestelman tolmlva satelllitti lahettaa ns. Ll-signaalia 
1575,42 MHz:n kantoaaltotaajuudella. TStS taajuutta merkitaan myos 
154fo. missa fo-10,23 MHz. LIsaksi satelllttit lahettavat l-2-signaalia 
1227,6 MHz:n kantoaaltotaajuudella. ell 120fo. Satelliitissa suoritetaan 

25 naiden signaalien modulointi ainakin yhdeUd valesatunnaissekvenssiila. 
Kullakin satelintilla tama valesatunnaissekvenssi on erilainen. Modu- 
loinnin tuloksena muodostuu koodimoduloitu laajakaistasignaall. 
Kaytetty modulointitekniikka mahdollistaa sen, etta vastaanottimessa 
pystytaan erottamaan eri satelliittien lahettamat signaalit, vaikka 

30 lahetyksessa kaytettavat kantoaaltotaajuudet ovat olennalsesti samat. 
Tasta modulointitekniikasta kaytetaan nimitysta koodijako- 
monikayttotekniikka (CDMA, Code Division Multiple Access). Kussakin 
satelliitissa LI -signaalin moduloinnissa kSytetS^n 

valesatunnaissekvensslna mm. ns. C/A-koodia (Coarse/Acquisition 

35 code), jona kaytetaan Gold-koodia. Jokainen GPS-satelliitti lahettaa 
signaalia kayttamalla ykslloilista C/A-koodia. Koodit muodostetaan 
kahden 1 023-bitttsen binSarisekvenssin mpdulo-2 summana. 
Ensimmainen bina^risekvenssi G1 on muodostettu polynomilla 
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ja toinen binaarisekvenssi G2 on muodostettu viivastamalla 
polynomia X^^+X^+X^+X^+X^+X^+l siten, ©tta kullakin satelliitilla viive 
on erilainen. Tama jarjestely mahdollistaa sen, atta eri C/A-koodit 
voidaan muodostaa samanlaisella koodlgeneraattorilla. C/A-koodit ovat 

5 siis binaarikoodeja, jolden kellotusnopeus (Chipping rate) GPS- 
jarjestelmassa on 1,023 MHz. C/A-koodi kasittaa 1023 alibittia (Chip), 
jolloin koodin toistoaika on 1 ms- L1-signaalin kantoaaltoa moduloidaan 
viela navigointi-informaatiolla 50 bit/s bittinopeudella, Navigointi- 
informaatio kSsittaa tietoa satelliitin "terveydentilasta" (health), radasta, 

10 aikatietoa jne. 

Satelliitit tarkkallevat laittelstonsa kuntoa toimintansa alkana. Satelliitit 
voivat kaytlaa esim, ns. vahtikolratoimintoja joidenkin laitteistoon 
mahdoHisesti tulleiden vikojen havaitsenniseen ja ilmoittamiseen, 

15 Virheet ja toimintahalrlot voivat olla hetkellisia tai pidempialkat$ia. 
Ten/eydentllatietojen pemsteella voidaan mahdollisesti osa yirhelsta 
kompensoida tai jattaa vjkantuneen satelliain lahettanria infonm 
kokonaan huomiQimatta. Lisaksi tilanteessa, jossa usecimman kuin 
neljan satelliitin signaall on vastaanotettavlssa, voidaan 

20 teiveydentilatletojen perusteeila painottaa eri sateiliiteista 
vastaanotettua informaatiota eri tavalla. Taildin epaiuotettavalta 
vaikuttavien satelliittien mahdollisesti alheuttamia virheita mittauksiin 
voidaan minimoida. 

25 Satelliittien signaalien havaitsemiseksi ja satelliittien tunnistamiseksi on 
vastaanottimen suoritettava tahdistus, jossa vastaanotin etsii kulloinkin 
kunkin satelliitin signaalin ja pyrkii tahdistumaan tahan signaaliin, jotta 
signaalin mukana l§hetettava data voidaan vastaanottaa Ja demodu- 
loida. 

30 

Tunnetun tekniikan mukaisissa vastaanottimissa tahan tahdistukseen 
kuluva aika riippuu mm, vastaanotetun signaalin voimakkuudesta. 
Tyypililsesti, mita heikompi vastaanotettava signaall on, sita pidempaan 
taylyy integroida jokalsessa etsintaavaruuden (korrelaatio/taajuus) 
35 alkiossa mahdollisen signaalin havaitsemiseksi. Tyypillisesti 
ulkokayttdon suunnitelluissa tunnetun tekniikan mukaisissa GPS- 
vastaanottimissa satelliittien signaaliin tahdistuminen vie muutamia 
kymmenia sekunteja tai muutamia minuutteja, mikall vastaanotettavan 
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signaalin voimakkuus on suhteellisen suuri. luokkaa -120- -130 dBm, 
Kuitenkin tilanteessa, jossa paikannus suoritettaisiin sisatiloissa tal sel- 
laisessa paikassa, jossa vastaanotettava signaali vaimonee eslmerkiksi 
rakennusten tai muiden maastoesteiden valkutuksesta, tahdistusaJka 

5 kasvaa huomattavasti. Sisatiloissa tyypiHinen sigriaaHvoimakkuus on n. :v" 
-150 dBm, iolloin tahdistusprosessiin kuluva aika sekventiaailisella^^^ 
haulla olisi tunteja. Tama aika on selvasti liian pitka jarkevan ^^^^ ^ >^ 
paikannuksen suorittamlseksi. Kuvassa 2a on esitetty r 
esimerkinomaisesti erasta satelliitista lahetettavaa CDMA- signaalia ja 

10 vastaavasti kuvassa 2b on esimerkinomaisesti esitetty tata lahetettya 
signaaiia sen saapuessa vastaanottimeen. Vastaanotettava signaali on 
snrtotiella vaimentunut huomattavasti ja siihen on summautunut 
kohlnaa. 

15 Sljalnninmaaritysvastaanottimen on suorltettava tahdistus mm. silloln, 
kun vastaanotin kytketaan paalle ja myos tilanteessa, jossa vastaanotin 
ei ole pitkaan aikaan pystynyt vastaanottamaan minkaan satelliitin 
signaalia. Mm. kannettavissa laltteissa tailainen tilanne voi syntya 
helposti, koska laite liikkuu ja iaitteen antenni ei aina ole optimaalisessa 

20 asennossa satelliltteililn nahden, mika heikentaa vastaanottimeen 
tulevan signaalin voimakkuutla. Myos kaupunkialuellla rakennukset 
valkuttavat vastaanotettavaan signaallin ja lisaksi voi syntya ns. ; {it 

monitie-etenemista, jossa IShetetty signaali saapuu vastaanottimeen eri 
kulkureitteja, esim. suoraan satelliitista (llne-of-sight) ja lisaksi 

25 rakennuksista heijastuneena, Tama monitie-eteneminen aiheuttaa sen, 
etta sama signaali vastaanotetaan useina eri vaitieisina signaaleina. 

Sijainninmaaritysjarjestelylle on kaksi p^aasiallista tehtavdS: 

30 1. vastaanottimen pseudo-etaisyyden laskenta eri GPS- 

sateiliitteihin, ja 

2. vastaanottimen sijalnnin maaritys, jossa kaytetaan 

laskettuja pseudo-etalsyyksia seka sateliiittlen sijaintitietoa^ 
Satelliittien kulloinenkin sijaintitieto voidaah laskea 
35 satellHteista vastaanotettujen ephemeris- ja 

aikakorjaustietojen perusteella. 
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Etaisyyksia satelliltteihin nimltetaan pseudo-et&isyyksiksi, koska aika ei 
vastaanottimessa ole tarkastl tiedossa. Talloin sijainnin ja ajan 
maaritysta toistetaan, kunnes on saavutettu riittava tarkkuus ajan ja 
sijainnin suhteen, Koska aikaa ei tiedeta absoluuttisen tarkastl, on 
5 paikka ja aika selvitettava linearisoimalla yhtaloryhma jokalsta uutta 
iteraatiota varten. 

PsGudo-etaisyyden laskenta voidaan suorittaa mittaamalla eri 
satelllittien signaalien keskinaiset, naennaiset kulkuviiveet. Sen jalkeen, 
10 kun vastaanotin on tahdistunut vastaanotettuun signaaliin, suorltetaan 
signaalissa lalietetyn infonnaation selvlttaminen. 

Lahes kaikki tunnetut GPS-vastaanottimet kayttavat kon^elaatlomene- 
telmia etaisyyksien laskentaan, Sijalnninmaaritysvastaanottimessa on 

1 5 tallennettu tai generoidaan paikallisesti eri satelllittien 
valesatunnaissekvenssit. VastaanotetuHe sigoaalijle suoritetaan muuntp 
valitaajuudelle (Down Conversion), minka jalkeen vastaanotin suorittaa 
vastaanotetun signaalin keitomisen tallennetulla 

valesatunnaissekvenssilla. Kertolaskuh tuloksena muodostunut signaali 

20 integroidaan tai alipSast6suodatetaan, jolloin tuloksena saadaan tieto 
siita, onko vastaanotetussa signaalissa ollut jonkin satelllitin latiettama 
slgnaalL Vastaanottimessa suoritettava kertolasku toistetaan siten, etta 
kullakin kerralla siirret§&n vastaanottimeen tallennetun 
valesatunnaissekvenssin vaihetta. Oikea vaihe paatellaan korrelaatiotu- 

25 loksesta edullisesti siten, etta korrelaatiotuloksen ollessa suurin, on oi- 
kea vaihe loytynyt. Talloin vastaanotin on oikein taiidistunut vastaan- 
otettuun signaaliin. 

Sen jalkeen, kun koodiin tahdistuminen on suoritettu, suoritetaan viela 
30 taajuuden liienosaatd ja vaihelukitus. Tama kdrrelaatiotulos ilmaisee 
myos GPS-signaalissa lahetetyn infonmaation, eli kyseessa on 
demoduloitu signaali- 

Edella mainittu tahdistus ja taajuudens^tdprosessi on toistettava kul- 
35 lekin sellaisen satellirtin signaalille, jota vastaanottimessa vastaanote- 
taan. Tahan prosessiin kuluu siis runsaasti aikaa erityisesti tilanteessa, 
jossa vastaanotettavat signaalit ovat lieikkoja. Taman prosessin no- 
peuttamlseksi kaytetaan joissaWn tunnetun tekniikan mukalslssa vas- 
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taanottimissa useampia korrelaattoreita, jolloin samanaikaisesti voidaan 
etsia useampia korrelaatiohuippuja. Pelkastaan korrelaattoreiden 
lukumaaraa kasvattamalla ei kaytannon sovelluksissa tahdistus- ja 
taajuudensaatoprosessia voida nopeuttaa kovlnkaan paljon, koska 
5 korrelaattoreiden lukumaaraa ei voi kasvattaa rajattomasti. 

Joissakin tunnetun tekniikan mukaislssa GPS-vastaanottimissa on kay- 
tetly FFT-tekniikkaa vastaanotetun GPS-signaalin Doppler-siirtyman 
maarittamlseksl tavallisten korrelaattoreiden yhteydessa, Naissa 
10 vastaanottimissa kaytetaan korrelointia vastaanotetun signaaiin 

kaistanleveyden kaventamiseksi 10 kHz-SO kHzilin. Tania^^^^ m: 
kapeakalstainen signaali anatysotelaan FFT-algoritmellla : < 
kantoaaltotaajuuden maarittamfseksi. 

1 5 Kansainv&lisessa patenttihakemuksessa WO 97/1 4057 on esitetty 
GPS-vastaanotin ja menetelma GPS-signaalien kasittelemisekisi. Tassa 
julkaisussa esitetty vastaanotin kasittaa paaasiassa kaksi eriiiista vas- 
taanotinta, joista ensimmainen vastaanotin on tarkoitettu kaytettavSksi 
tilanteessa, jossa vastaanotetun signaaiin voimakkuus on rlittavan suuri 

20 ja toista vastaanotinta on tarkoitettu kaytettavaksi tjlanteessa, jossa 
vastaanotetun signaaiin voinnakkuus ei riita rlittavan tarkkaan sijainnin 
maaritykseen kaytettaessa enstmmaista vastaanotinta. Tassa toisessa 
vastaanottimessa suoritetaan vastaanotetun signaaiin digitointi ja tal- 
iennus muistivailneisiin, jolloin mydliemmin naita tallennetttuja signaalei- 

25 ta kasiteiiaan digitaalisessa signaalinkasittely-ykslkdssa. Digltaalinen 
signaalinkaslttely-yksikko suorittaa konvoluiitio-pperaatioita vastaanpte- 
tulie digltoidulle signaalllle. Ndiden konvoJuutio-operaatioiden tavoittee^ 
na on pseudo-etaisyyksien laskenta. Muistivalinelslln tallennetaan 
tyyplllisesti sadasta tuhanteen koodijaksoa (PM Frames), mikS vastaa 

30 100ms-1s mittaista signaalla. Taman jalkeen tutklttavan satelliitin 
koodia vastaava tallennettu koodi noudetaan vastaanottimen muistista 
kaytettavaksi vastaanotetun signaaiin analysolntiin. 

Vastaanottimessa poistetaan myos Doppler-siirtyma. Taman Doppler- 
35 siirtyman maara on selvitetty joko ensimmaisella vastaanottimella tai 
GPS-jarjestelmaan kuuluvasta tukiasemasta vastaanotetun tiedon pe- 
njsteella. Taman jalkeen suoritetaan perakkaisten kehysten koherentti 
summaus. Talle summauksen tuioksena saadulle datajoukolle suorite- 
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taan nopea Fourier-muunnos. Fourier-muunnostulokselle suoritetaan 
kertolasku, jossa toisena kertojana kaytetaan muistivalineisiin tallenne- 
tun referenssisignaalin Fourler-muunnoksen kompleksikoniukaattia. 
Talle kertolaskun tulokselle suoritetaan edelleen kaanteinen Fourler- 
5 muunnos, jolloin saadaan joukko korrelaatiotuloksia. Tassa julkalsussa 
on siis korrelaatio korvattu Fourier-muunnoksella, jolloin laskutoimitus- 
ten mSaraa on saatu pi©nennettya. Julkaisun mukaan menetelma no- 
peuttaa 1 0-1 00 -kertaisesti sijalnninmaaritysta verrattuna mainitun jul- 
kaisun hakemlshetkellatunnettulhin ratkaisuihln. 

10 

Nyt esilla olevan keksinnon erasna tarkoituksena on aikaansaada vas- ^ yi 
taanotin, jossa tahdistus lahetettyyn signaaliin voidaan suorlttaa nnerkit^ 
tavassa maSrin nopeammin kuin tunnetun teknilkan mukaisissa vas- 
taanottimlssa myos heikoilla signaalellla. Keksint6 soveltuu erityisesti 

15 sijainninmaaritysvastaanottimissa kaytettavaksi, mutta myos muissa 
vastaanottimissa, edullisesti CDMA-vastaanottimissa, jolssa vastaanot- 
timen on tahdistuttava hajaspektrisignaaliin. Keksinto pemstuu sHhen 
ajatukseen, etta hyodynnetaan konvoluutlon ja aika- 
taajuusmijunnoksen, kuten nopean Fourier-muunnoksen, 

20 taajuusanalyysiominaisuuksia tehokkaammin, jolloin samanaikalsesti 
voidaan tutkia suurempi joukko informaatiota kuin tunnetun tekniikan 
mukaisissa Fourier-muunnosta kayttavissa ratkaisuissa, Nyt esilla ole- 
van keksinnon mukaiselle menetelmalle on tunnusomaista se, mlta on 
esitetly oheisen patenttivaatimuksen 1 tunnusmerkklosassa* Nyt esilla 

25 olevan keksinnon mukaiselle vastaanottimelle on tunnusomaista se, 

mlta on esltetty oheisen patenttivaatimuksen 17 tunnusmerkkiosassa. • • 
Nyt esilla olevan keksinnon mukaiselle elektroniikkalaitteelle on vieia 
tunnusomaista se, mlta on esitetty oheisen patenttivaatimuksen 32 
tunnusmerkkiosassa. 

30 

Nyt esllia olevalla keksinnolla saavutetaan merkittavia etuja tunnetun 
tekniikan mukaislin menetelmlin ja vastaanottimiin verrattuna. Keksin- 
non mukalsella menetelmalla voidaan erityisesti sijainninmaaritysvas- 
taanottimissa nopeuttaa vastaanottimen tahdistusta vastaanotettavaan 
35 signaaliin. Talldin sijainnin maaritys voidaan suorittaa huomattavasti 
nopeammin my6s tilanteissa, jolssa vastaanotettava signaali on heikko, 
eslmerkiksl sisatlioissa tai kaupunkialueilla, jolssa slgnaalin etenemis- 
tlella on signaalia merkittavasti vaimentavia esteita. Keksinnon mukai- 
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sella menetelmalla voidaan saavuttaa merkittava parannus 
tahdistusnopeuteen verrattuna tunnetun tekniikan mukaisiin menetel- 
miin. Keksinnon mukainen vastaanotin voidaan toteuttaa suhteellisen 
pienella komponenttimaaralla seka kokonaisenergiankulutus voidaan 
5 pitaa kohtuullisena, jolloin keksinto sovettuLi erityisen hyvin 
kaytettavaksi kannettavissa laitieissa. TSIIoln 

sijainninmaaritysvastaanotin voidaan toteuttaa myos langattoman 
viestimen yhteyteen. 

10 Nyt esilia olevaa keksintda selostetaan seuraavassa tarkemmin viitaten 
samalla oheisiin piirustuksiin, joissa 



15 



kuva 1 eslttaa pelkistettyna lohkokaaviona vaistaanotinta, jossa 
keksinnon mukaista menetelmaa voidaan soveltaa, 

kuva 2a esittaa yksinkertaistettuna satelliitista lahetettavaa signaaiia, 



20 



kuva 2b esittaa esimerkkiS kuvan 2a yksinkertaistetusta signaalista 
vastaanottimessa, 

kuva 3 esittaa keksinnon eraan eduHisen suorltusmuodon 
mukaisen menetelman korreiolntlvailietta pelkistettyna 
kaaviona, 



25 kuva 4 



esittaa keksinnon eraan edullisen suoritusmuodon 
mukaisen menetelman analyysivaihetta pelkistettyna 
kaaviona, 



30 



kuva 5 esittaa keksinnon eraan edullisen suoritusmuodon 
mukaisen menetelman summausvaihetta pelkistetysti, 



kuva 6 esittaa keksinnon eraan edullisen suoritusmuodon 
mukaisen menetelman paattelyvaihetta pelkistetysti, 



35 kuva 7 



esittaa keksinnon eraan edullisen suoritusmuodon mukaista 
eiektroniikkalaitetta pelkistettya lohkokaaviona. ja 
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kuvaS esittaa erasta tunnetun tekniikan mukaista korrelaattoria, 
jota voidaan kayttaa keksinnon eraan edullisen 
suoritusmuodon mukalsen menetelman yhteydessa. 

5 Kuvan 1 vastaanottimessa 1 vastaanotettava signaali muunnetaan so- 
pivimmin valitaajuudelle muunninlohkossa 2. Tassa vaiheessa signaali 
kasittaa sinansS tunnetusti kaksi komponenttla: I- ja Q-komponentit, 
joiden valilla on 90° vaihe-ero, Nama valitaajuudelle muunnetut analo- 
glset signaalikomponentit digitoidaan digitointilohkossa 3 ja johdetaan 

10 kertojalohkoon 4. Kertojalohkossa 4 digitoidun signaalln I- ja Q-kompo- 
nentit kerrotaan numeerisesti ohjatun oskillaatlorin 5 (NCO, Numerically 
Controlled Oscillator) muodostamalla signaalllla, Tama numeerisesti 
ohjatun osklllaattorin 5 signaali on tarkortettu korjaamaan Doppler-siir- 
tymasta ja vastaanottimen 1 paikallisoskillaattorin (ei esitetty) taajuus- 

15 virheesta Johtuva taajuuspoikkeama. Tata numeerisesti ohjattua oskil- 
laattoria 5 kaytetaan myos keksinnon mukaisen menetelman yhtey- 
dessa, kuten myohemmin tassa selityksessa esitetaan/Kertojato^^ 
muodostama signaali johdetaan tahdistuslohkoon 6, jossa kekslnnGn 
mukaista meneteimaa sovelletaan. T§ma tahdistuslohko 6 pyrkii siis : 

20 loytamaan satelliitin koodivaiheen ja taajuuspolkkeaman kaytettavaksi 
tahdistuksen jalkeislssa toiminnoissa. Myos tata kuvataan jaljempana 
tassa selityksessa. Ohjauslohkolla 7 ohjataan tahdistuksen aikana 
skannauslohkoa 8, jonka avulla numeerisesti ohjatun osklllaattorin 5 
taajuutta saadetaan tan/lttaessa. Ohjauslohko 7 ohjaa ensimmaista 

25 kytkinta 9 kytkemaan joko taman skannauslohkon 8 muodostaman 
signaalln tahdistuksen aikana numeerisesti ohjatulie oskillaattoriile 5 tai 
koodin selvityslohkon 11 muodostaman ohjaussignaalin numeerisesti 
ohjatulie oskillaattoriile 5 sen jalkeen, kun tahdistus on saavutettu. 
Toiseila kytkimella 10 ohjataan tata seurantalohkon 11 toimlntaa, Tama 

30 seurantalohko 11 muodostaa osan sinansa tunnetusta 
koodlvaihelukitusta silmukasta ja kantoaaltovaiheiukltusta silmukasta 
(ei esitetty). 

Kayttojannitteiden kytkemlsen jalkeen tai tilanteessa, ibssa vastaanotin 
35 1 ei ole pystynyt vastaanottamaan GPS-satelliittlen signaalia pitkaan 
aikaan, suoritetaan vastaanottimessa 1 kaksidimensioinen etsintavaihe 
kullekin satelliitille, joiden signaalia vastaanotetaan. Tassa 
kaksidimensioisessa etsintavaiheessa tarkoituksena on seivittaa kunkin 



09/1 2 '99 18:26 FAX 03 2886262 



TAMPEREEN PATENT 12)010 



satelliitin kantoaaltotaajuus seka koodivaihe. Tahan^^^^^^. v. 
kantoaaltotaajuuteen vaikuttaa siis satelliitin lilkkumisesta johtuva; 
Doppler-siirtyma seka vastaanottimen paikallisoskillaattorin 
epatarkkuudet. Taajuusepavarmuus vol olla varsin suurl, jopa ± 6 kHz. 

5 jolloin vastaanottimen 1 on suoritettava etsinta n. 12kHz:n 
taajuusalueelta varsinaisen lahetystaajuuden (L1 = 1575,42 MHz) 
suhteen. Vastaanotin 1 el myoskaan tieda tarkkaa koodivaihetta, jolloin 
vastaanottimen on suoritettava my6s koodivaiheen selvitys 1023 
mahdollisesta koodlvalheesta. Tail5in saadaan kaksidimensloinen 

10 etsintaprosessi, Jossa yksi dimenslo on taajuuspoikkeama 12 kHz:n 
alueella ja toinen dimensio on koodivaihe 1023:sta eri|alsesta koodi- 
vaiheesta. Keksinnon er^an edulllsen suoritusmuodon mukaisessa 
menetelmassa voidaan kerrallaan tutkia n. 500 Hz:n taajuusalue, jolloin 
menetelmaa toistetaan tarvitlaessa 24 kertaa koko tutklttavan 12 kHz:n 

15 taajuusalueen kattamiseksi. On selvaa, ett§ tassa selltyksessa 
kaytettavat esimerkkiarvot ovat vain k^ksintd^ selyentavina, mutla eivSt 
rajolttavina esimerkkeina. KeksintOa voidaan soveltaa mutssakin kuin 
GPS-jarjestelmlssa, jolloin mm. mainitut taajuusarvot. koodivaiheet 
sek& koodien lukumMra vol valhdella. 

20 

Selostetaan seuraavaksi keksinnon eraan edullisen suoritusmuodon 
mukaisen menetelm^n toimintaa kuvan 1 mukaisessa vastaanotti- 
messa 1 . Tahdistuksen kaynnistamiseksi skannauslohko 8 asettaa nu- 
meerisesti olijatun oskillaattorin 5 taajuuden siten, etta vastaanottimella 

25 1 vastaanotetaan edullisesti taajuusalueen pienimpia taajuuksia, tassa 
esimerkissa 1575.414 MHz — 1575,4145 MHz. Vastaanotin voi 
maarittaa aloitustaajuuden myos siten, ettS vastaanotin kayttaa hyvaksi 
esim. aikaisemmin seMtettya sijaintitietoa ja/tai almanakkatietoa, jolloin 
sijainninmaaritysta voidaan vieia nopeuttaa. Vastaanotettavasta sig^ 

30 naalista tallennetaan naytteita naytevektoreiden muodostuslohkossa 12 
edullisesti komplekslsiksi naytevektoreiksi Pk(i}. Pk(2)...pk(N), joissa 
kussakin on 1023 naytetta tassa edullisessa suoritusmuodossa. 
Naytevektoreiden muodostuslohkon 12 n^ytteiden tallennustaajuus 
tassa edullisessa suoritusmuodossa on olennEusesti sama kuin 

35 aliblttlen kellotusnopeus, eli n. 1 023 000 naytetta sekunnissa. 
Naytevektorit ovat jatkuvia siten, etta seuraava naytevektori jatkuu 
ajallisesti edellisen naytevektorin jalkeen, eli naytevektorin viimersen 
naytteen ja seuraavan naytevektorin ensimmaisen naytteen aikaero on 
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olennaisesti sama kuin naytevektorin perakkaisten naytteiden valinen 
aikaero. Nama 1023 naytetta vastaavat siis n. 1 ms:n mittalsta 
signaalia, jolloin aika-taajuusmuunnoksessa taajuusalue on n. 1 kHz, 
Josta osa voidaan hyodyntaa. Naytevektoreklen muodostu3vaihetta on 
5 esitetty viitteella 101 oheisessa kuvassaS- 

Naytevekloreita on edullisesti N kappaletta, missa N on sopivimmin 
kahden potenssK Lisaksi naytevektoreiden muodostus toistetaan 
keksinnon edullisessa suoritusmuodossa K kertaa, kuten myShennmIn 

10 tassa selityksessa ositetaan. Alaindeksilla k ilmaistaan seuraavissa 
merkinnolssa eri toistokertoja. Naytevektoreiden pk(1), pk(2).»Pk(N) 
lukumaaran N arvoa maaritettaess& tulee huomioida GPS-jarjes- 
telmassa se, etta signaaliin on moduloitu informaatiota bittinopeudella 
50 bittia/s binaarisena vaihemodulaatlona. Tolnen tata naytevekto- 

15 reiden Pk(1), pic(2)...Pk(N) lukumaaraa N rajoittava teklja on 
vastaanottimen paikallisoskillaattorin taajuusst^^^^ 

Naytevektoreiden mucdostusvaiheen llsSksi keksinnon mukaisessa 
tahdistusmenetelmassa suoritetaan korrelointlvaihe 

20 korrelaatiofunktiomatriisin muodostamiseksi. 

Tama korrelointlvaihe voidaan suorlttaa osittain jo naytteenoton aikana, 
tai sen iaikeen kun N kappaletta naytevekloreita Pk(1), Pk(2)-..Pk(N) on 
muodostettu. Jos korrelointlvaihe suoritetaan esim. siten, etta kunkin 

25 naytevektorin tallennuksen jaikeen lasketaan sille aika- 
taajuusmuunnos. kuten nopea Fourier-muunnos (PfT), voidaan samaa 
aika-taajuusmuunnlnta kayttaa kaikissa N kappaleessa naytevektoreita 
Pk(1). Pk(2)...pK(N). Jos sen sijaan korrelointlvaihe suoritetaan N nayte- 
vektorin tallennuksen jaikeen. on kaytettava joko kullekin naytevektoriile 

30 omaa aika-taajuusmuunninta, tai aika-taajuusmuunnokset suoritetaan 

eri naytevektoreille perakkain samassa aika-taaiuusmuuntimessa. 

Kuvassa 3 on esitetty menetelman korrelointivalhetta, jossa naytevekto- 
reista Pk{1), Pk(2)...pk(N) muodostetaan korrelaatiofunktiomatriisl Qx k . 

Vaikka jatkossa tassa selityksessa paaasiassa kaytetaan esimerkldna 
35 aika-taajuusmuunnoksesta Fourier-muunnosta ja kaanteisesta 
muunnoksesta, eli taajuus-aikamuunnoksesta kaanteista Fourier- 
muunnosta. on selvaa, etta nyt esilla olevaa kekslntoa ei ole rajoitettu 
ainoastaan naihin eslmerkkeihin. 
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Kullekrn naytevektorille Pk(1). pk(2)...pk(N) suoritetaan diskreetti Fourier- 
muunnos 102, sopivimmin nopea Fourier-muunnos, eli FFT-muunnos: 

5 Pk(i) = FFT(p^(i)), missai=1,..., N (1) 

Tata esittavat lohkot FFT1, FFT2 FFTN kuvassa 3. Kaytannon las- 

kutoimituksissa kSytetaan edulllsesti 1024:aa arvoa, koska tallQin 
diskreetti Fourier-muunnos on kaytanndn sovelluksissa huomattavasti 
10 tehokkaammin toteutettavissa (FFT-a!goritmllla) kuin kaytettaessa 
1023:a afvoa. Eras tapa taman aikaansaamiseksl on lisata 
ylimaarainen nolla 1024:ksi alkioksi. Talla on vahainen vaikutus 
muunnostulokseen . 

15 Vastaanottimessa on tallennettuna sopivimmin kaikkien GPS-jarjestel- 
man sateiiiittien (ei esitetty) C/A-koodia vastaavat vertailukoodit r{x), 
missa X virttaa satelilitin tunnukseen, eslm. x on vSlillg^ 1 — 36. 
Vertaiiukoodeja el valttamatta tarvltse tallentaa, vaan ne voidaan myos 
generoida vastaanottimessa. Korrelaatiovaiheessa valitaan tal 

20 generoidaan kulioinkin sen satelilitin vertailukoodi, jonka laliettamaan 
signaaliin vastaanotinta tahdistetaan. Vertailukoodi on kflSnnetty 
ajallisesti takaperin. Talle takaperoiselle yertailukoodille. jota kuvassa 3 
on merkitty viitteella f(x) . suoritetaan diskreetti Fourier-muunnos 103, 

sopivimmin nopea Fourier-muunnos, eli FFT-muunnos: 

25 

R(x) = FFT(r(x)) (2) 

Takaperoinen vertailukoodi r(x) ja/tai sen FFT-muunnos on voitu 

tallentaa jo etukateen vastaanottimen muistivaiineisiin, tai se 
30 muodostetaan vertailukoodista r(x) tahdistuksen yhteydessa. 

Seuraavaksi korrelaatiovaiheessa suoritetaan kunkin naytevektorin Pk(i) 
Fourier-muunnostuloksen P^O) ja tai<aperoisen vertailukoodin f(x) 

Fourier-muunnoksen R(x) valinen kertolasku 104: 

35 _ 

M,,,(i)=R(x)-Pk(i) (3) 
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Naille kertolaskujen tuloksille suoritetaan viela kaanteinen Fourler- 
muunnos 105, jolloin tuloksena saadaan vertailukoodin r(x) ja 
vastaanotetun signaalin ristikorrelaatio kaikilla mahdoHisllla 
5 kokonaislukuviipeilla (1023kpl). 

m,,,(i) = iFFT(M,,,(i)) (4) 

Tama perustuu siihen. etta aikatason signaalien konvoluution Fourier- 
10 muunnos vastaa Fourier-muunnettujen, eli taajuustasoon muunnettujen 
aikatason signaalien, kertolaskua. Kun llsaksl kaytetaan takaperoista 
vertailukoodia, voidaan Fourier-muunnoksella suorittaa nopea 
diskreettiaikainen korrelaatio. Ristikorrelaatiotulos kasittaa tassa edul- 
lisessa esimerkissa talloin 1023 alkiota. NAIsta eri nSytevektoreista Pk(l) 
15 muodostetuista ristikorrelaatlotulokslsta mx,k(i) muodostetaan korrelaa- 
tiofunktlomatrilsi C^^^, jossa rivien maara on naytevektoreiden 

lukumaara N. 

On selvaa, etta vertailukoodin ajaliisen kgidntamisen sijasta voidaan 
20 naytevektoreista Pk(i) muodostaa takaperoisia naytevekloreita pk(i), 

joHoin edelia esitetyissa laskutoimKuksissa kaytetaan vertailukoodia r(x) 
suoraan ja takaperoisia naytevektoreita. EraassSi edullisessa 
suoritusmuodossa ei tarvitse tehda kumpaakaan edella mainituista 
kaantamisista, eli voidaan kaytt^a vertailukoodia r(x) ja naytevektoreita 
25 pk(i) seilaJsenaan. Tama perustuu korrelaattoteoreeman ominaisuuden 
hyddynt^iseen, joka ilmaisee, etta kahden aikadiskreetin funktion z^, 
Za valinen ristikorrelaatio con'(zi, Za) voidaan muodostaa taajuustasoon 
muunnettujen funktioiden taajuus-aikamuunnoksen avulla. Tama 
voidaan esittaa seuraavalla kaavalla: 

30 

corr(z,.Z2)oc|FFT(z;.Z2) (5) 

Talloin keksinnon eraassa edullisessa suoritusmuodossa voidaan 
naytevektoreiden Pk(i) ja vertailukoodin r(x) valinen ristlkorrelaatfo 

35 muodostaa esim. suorittamalla naytevektoreille Pk(i) ja vertailukoodrlle 

r(x) Fourier-muunnos (soplvlmmln nopea Fourier-muunnos FFT) 
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naytevektorin Pk(i) Fourier-muunnoksen Pk(i) ja vertailukoodin r(x) 
Fourier-muunnoksen R(x) muodostamlseksi. muodostamalla 
naytevektorin Fourier-muunnoksesta Pk(i) kompleksikonjugaatti P*k(i) . 
kertomalla muodostettu kompleksikonjugaatti P\{\) naytevektorin 
Fourier-muunnQksellaR(x) , Ja suorittamalla kertolaskutulokselle 
kaanteinen Fourier-muunnos. Vaihtoehtoisesti voidaan muodostaa 
vertailukoodin Fourier-muunnoksesta R(x) kompleksikonjugaatti R (x), 
jolloln se kerrotaan naytevektorin Fourier-muunnoksella minka 
Jalkeen suoritetaan kertolaskutulokselle kSanteinen Fourier-muunnos. 



Tassa yhteydessa on syyta korostaa sita. etta edelllsessa kappaleessa 
esttetty menetelmS naytevektorelden Pk(i) ja vertailukoodin r(x) valisen 

ristlkorrelaation laskemiseksl johtuu korrelaatlon ja konvoluution 
pemsominaisuuksista Ja niiden vaiisesta lahelsesta riippuvuudesta, 

15 Jolloin funktion kaantamisella takaperoiseksi aikatasossa on 
kaytannOltlsestI katsoen matemaattinen riippuvuus 

kompleksikonjugaatin muodostamiseen taajuustasossa. Tata asiaa on 
tarkemmin kasitelty mm. julkaisussa: "Digital Signal Processing - A 
Practical Approach", Emmanuel C. Ifeachor and Barrie W. Jervis, 

20 Addlson-Wesley Publishing Company 1993, ISBN 0-201 -5441 3-X. 
kappale 4: "Correlation and Convolution", johon tassa yhteydessa 
viitataan. On syyta vlelS mainita, etta kokonalsuudessaan taman 
kekslnnon soveltamisen kannalta ei singnsa ole merkitysta silia, mita 
menetelmaa kaytetaan ristikon-elaatiotuloksen aikaansaamiseksi. 



Korrelaatiovaiheessa muodostetun korrelaattofunktiqmatriisin Cx.k nyit 
eslttavat vastaanotetun signaalin ja verfcailukoodin ristikorrelaatiota eri 
vaihe-eroilla yhden mlllisekunnin valein otettuna. Kaavana 
korrelaatiofunktiomatriisi voidaan talldin esittaa seuraavasti: 



m 



x,k 



0) 

m^(2) 
mx.k(N) 



(6) 
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Seuraavassa vaiheessa eli analyysivaiheessa kaytetaan 
korrelaatiofunktiomatriisin C^.k transpoosia 106, jossa rivit esittavat 

aikatasossa signaalin naytteita siten kuin tunnetun tekniikan 
mukaisessa korrelaattorissa. Kukin rivi vastaa tiettya koodivaihe-eroa 
vastaanotetun signaalin ja vertailukoodin valilld. Taman 
korrelaatjofunktiomatiilsin Cx.k transpoosin kullekin rlville suoritetaan 
Fourier-muunnos 107 koherentin etsintamatrjisin Ax,i< 
muodostamiseksi, jolloin voidaan suorittaa taajuusanalyysi todeilisen 
taajuussiirtyman selvittamiseksl. Tata eslttaa olieinen kuva 4, 



Kaytannon sovellukslssa ei korrelaatiofunktiomatriisista tarvitse 
erikseen muodostaa transponoitua matriisia, vaan tallennetun 
korrelaatiofunktiomatriisin Cx,k alkiot luetaan muistista 16 (kuva 7) eri 

suunnassa, eduliisesti sarakkeittain. 

Korrelaatlofunktiomatriisi Cx,k voidaan muodostaa myos eslm. sinansa 
tunnetusti sovitettuja suodattimia (matched filter) kayttamalla. 
Esimerkiksi patentissa US-5,373,531 on esitetty eras sovitetuilla 
suodattimiila toteutettu korrelaattori 25 (kuva 8), jota voidaan soveltaa 
korrelaatiofunktiomatriisin Cx,k muodostamiseksi. Sovitettu suodatin on 

muodostettu kolmesta siirtorekisterista siten. etta ensimmaiseen 
siirtorekisteriiri 26 jolndetaan vastaanotetun signaalin l-komponentti, 
toiseen siirtorekisteriin 27 jolidetaan vastaanotetun signaalin Q- 
komponentti. ja kolmanteen siirtorekisteriin 28 joiidetaan satelliitin 
koodia vastaava vertailukoodl. Siirtorekistereiden 26, 27, 28 pituus on 
sopivimmin naytevektorin pituus, eli 1023 bittia GPS-jarjestelmassa. 
Kun yksi naytevektori on siirretty siirtorekisteriin 26. 27, suoritetaan 
korrelointi seuraavasti. l-komponentin sisaltavan siirtorekisterin 26 
bittien ja vertailukoodin sisaltavan siirtorekisterin 28 bittien vaiilia 
tehdaan ehdoton-ei-tai -operaatio (Exclusive NOR) ensimmaisessa 
vertailulohkossa 29 ja Q-komponentin sisaltavan siirtorekisterin 27 
bittien ja vertailukoodin sisaltavan siirtorekisterin 28 bittien valilia 
tehdaan ehdoton-ei-tai -opersiatio (Exclusive NOR) toisessa 
vertailulohkossa 30. Operaation tuloksena saadaan bittikohtainen 
korrelaatiotieto: jos vastaanotetun signaalin (I- ja Q-kompbnenteittaIn 
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tarkasteltuna) bitti ja vastaava vertailukoodin bitti tasmaavat. tuloksena 
on blna^riarvo 1 , muussa tapauksessa tuloksena on binaariatvo 0. Jos 
signaall korreloi taysin vertailukoodin kanssa. kaikki bitit ovat arvossa 1 . 
Jos mikaan bitti ei korreloi, kaikki bItIt ovat arvossa 0. Korrelaation 

5 tuloksena saadut arvot summataan ensimmalsessa summauslohkossa 
31 (l-arvot) ja toisessa summauslohkossa 32 (Q-arvot), Jolloin 
tuloksena saadaan korrelaatioarvo mx,k(i) tietylla koodiviiveella tiettyna 
alkana. Yhden vertaiiun jaikeen suoritetaan vertailukoodin siirto yhdelia, 
minka jaikeen muodostetaan uusi korrelaatiotulos. Kukin edella 

10 esitetylla menetelmalia muodostettu korrelaatiotulos vastaa yhta 
korrelaatiofunktiomatriisin Cx,k "vin alkiota (yhteensa 1023). Edella 

mainittuja vaiheita toistetaan tarvitlaessa kullekin nayteveklorille, jolloin 
kullakin laskentakerralla saadaan korrelaatiofunktiomatriisin Cx,k yksi 

rivi. 

15 

Kuten edella kuitenkln todettiin, GPS-jarjestelmassa signaalia 
moduloidaan 50 bit/s signaalllla. mika rajoittaa lukumaaran N arvoa 
kaytannon sovelluksissa. TallOin lukumaara N on valittava sopivlmmin 
siten, etta modulointi ei olennaisesti valkuta analyysiin. Lisaksi optimi 

20 N-arvoon vaikuttaa se, mita ikkunafunktiota Fourier-muunnoksessa 
kaytetaan. Valitsemalla lukuarvoksl esim. N = 32 saadaan 
vastaanottimen kohinakaistan leveydeksl luokkaa 30 Hz, joka on viela 
hieman liian suuri sellaisten signaalien Hmaiaemiseen, jolden 
voimakkuus on luokkaa -150 dBm vastaanottimessa. Taman vuoksi 

25 tahdistuslohkossa 6 suoritetaan vieia optkinaalinen epakoherentti sum- 
mausvaihe, jossa signaaJIkohinasuhdetta parannetaan. 

Epakoherentin summausvaiheen toteuttamiseksi toistetaan 108 edella 
esrtettyja naytevektoreiden muodostusvaihetta, korrelointivaihetta ja 

30 analyysivaihetta K kertaa (kuva 5). Tama toistojen lukumaara K valitaan 
eduHisesti siten, etta slgnaall-kohinasuhdetta saadaan riittavasti 
parannettua, kuitenkin kohtuullisessa ajassa. Kullakin analyysivaiheen 
suorituskerralla muodostetaan yksi koherentti etsintamatriisiAx,k> jo"© 
suoritetaan epakoherentti summaus epakoherentin etsintamatrlisin Sx 

35 muodostamiseksi. Epakoherentti etsintamatriisi muodostetaan 
edullisesti seuraavasti. Kunkin koherentin etsintamatriisin Ax,k 
kompleksisista alkioista ax^K{i,j) lasketaan edullisesti suuruus (itseis- 
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arvo) tai muu lukuarvo, kuten alkion suuruusarvon toiseen potenssiin 
korotus. Kunkin epakoherentin etsintamatriisin vastinalkioista lasketut 
lukuarvot summataan, ell suoritetaan matriisien yhteenlasku, joka 
voidaan esittaa kaavalla: 



K 

k=1 



Ja^,(l023,1)| ... |a^K(l023,N)|_ 



(8) 



Epakoherentti etsintamatriisi voidaan kaytannOn sovelluksissa 
muodostaa ainakin kahdella tavalla. Ensimmaisessa 

10 toteutusvaihtoeJidossa tallennetaan kullakin toistokerralla muodostettu 
koherentli etsintamatriisi. Tarvittavien toistokertojen jalkeen suoritetaan 
epakoherentti etsintamatriisi summaamalla vastinalkiot kaavan 8 
mukalsesti. Tassa toteutusvaihtoeJidossa tarvitaan muistia kalkklen 
koherenttien etsintamatrllsien alkioiden tallentamiseen. Toisen 

1 5 toteutusvaihtoehdon mukaan lasketaan ensin yksi koherentti 
etsintamatriisi, jonka arvot kopioidaan summamatriisin alkioiksi. 
Kullakin seuraavalla toistokerralla muodostetaan koherentti 
etsintamatriisi. jonka arvot summataan ep&koherentin etsintamatriisin 
vastlnalkioiden kanssa. TassS valhtoehdossa suoritetaan 

20 vastinaJkiolden summaus siis joka toistokerralla. Talloln tallennetaan 
vain yksi koherentti etsintamatriisi, jolloln mulsMa tarvitaan vahemmSn 
kuin ensimmaisesa valhtoehdossa. 

Sen jalkeen, kun tarvittavat toistokerrat on suoritettu, selvltysvaiheessa 
25 tutkitaan taman epakoherentin etsintamatriisin alkioiden Sx(U) 

arvoja ja pyrltaan loytamaan 110 sellalnen arvo, joka ylittaa ennalta 
asetetun kynnysarvon ja on selvasti multa arvoja suurempi. Jos 
tallainen arvo loytyy, ilmaisee se koodlvaihe-eron seka taajuus- 
polkkeaman, koska Kyseessa on todennakoisesti satelliitin lahettama 

30 signaali. Mikali signaali ei ole satelliitin lahettamaa, vaan kohinaa tai 
muuta satunnaista hairiOsignaalia, ei merkittavia konrelaatiohulppuja 
pitaisi muodostua. Koodivaihe-ero ilmene© taman suurimman arvon 
rivi-indeksista ja vastaavasti taajuuspoikkeama llmenee sarake- 
indeksista. Kuvan 6 eslmerklssa on esltetty eras tilanne. jossa selvasti 

35 multa arvoja suurempi arvo on loytynyt. Jos sen sijaan epakoherentista 
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etsintamatriisista ei loydy tallaista arvoa, eli tutkitulla taajuusalueella 
ei todennakoisesti vastaanotettu haettavan satelliitin lahettamaa 
signaalia, muutetaan tutkittavaa taajuuskaistaa ja suoritetaan edGlla 
esitetyt vaiheet epakoherentin etsintamatriisin muodostamiseksi. Talla 
5 menetelmaila saadaan koko tutkittava 6 kHz:n alue kaytya lapi 
toistamalla edella esitetyt vaiheet tarvittavan monta kertaa. 

Tarvittaessa voidaan edella esitetyt vaiheet toistaa koko tutkittavalle 
taajuuskaistalle ja tallentaa eri toistokerroilla muodostetut epakohererrtit 
10 etsintamatriisit tai vain mahdolliset huippukohdat, ennen kuin etsitaan 
suurinta korrelaatiohuippua, Talloin voidaan pienentaa virhetuikintojen 
mahdollisuutta mm. sellaisessa tilanteessa. jossa kynnysarvo on 
asetettu iiian pieneksi ja hairidsignaali vol aiheuttaa virhetulkinnan. 

15 Jos olisi saatavilla laittelta, jotka pystyisivat suorittamaan koodivalheen 
ja taajuuspoikkeman signaaleilla, joiden voimakkuus on luokkaa 
-150 dBm, olisi suoritettava taajuusalueen skannaus kymmenlen 
hertzien levyisissa osissa, jotta vaihetta ja taajuuspoikkeamaa voitaisiin 
edes yrittaa selvtttaa. Nyt esllla olevalla kekslnnolla voidaan suorittaa 

20 taajuusalueen skannaus satojen hertzien levyislssa osissa, sils 
huomattavasti nopeammin kuin tunnetun tekniikan mukaisilla 
vastaanottimilla. 

Sen jalkeen kun oikea taajuuspoikkeama ja koodivaihe on seM'tetty, 
25 voidaan vastaanotin asettaa seurantatllaan (tracking)- Helkoimmilla 

signaaleilla ei datavastaanotto onnistu, vaan on sin^nsa tunnetusti 

turvauduttava esim. matkaviestinverkon kautta saatuun dataan. 

Etaisyysmittaus on edelleen mahdollista alentuneella tarkkuudella, 

Seurantatila asetetaan vaihtamalla kytkimet (kuva 1) toiseen asentoon, 
30 jolloin vastaanotettu inforniaatlo johdetaan seurantalohkoon 1 1 , jossa 

muodostetaan myos takaislnkytkenta numeerisesti ohjatun oskillaattorin 

5 taajuuden hienosaatoa varten. 

Sijainnin laskemlseksi vastaanotin suorittaa signaalin vastaanoton so- 
35 pivimmin vahintaan neijasta satelliitista vastaanotetun signaalin perus- 
teella. Talloin edella esitetty tahdistus toistetaan tarvittaessa kunkin 
satelliitin signaalia varten, jolloin vertailusekvenssiksi r(x) valitaan kul- 
loinkin sen satelliitin koodi, johon tahdistus suoritetaan. 
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Suuri osa menetelman toteuttamisessa tarvittavista lohkoista voidaan 
toteuttaa esim. digitaalisessa signaalinkasittely-yksikossa (Digital 
Signal Processor, et esitetty). FFT-muunnosten suorittamisessa voi- 
6 daan kayttaa joko laittelstopohjaisia ratkaisuja tai digitaalisen signaalin- 
kasittely-yksikon ohielmalllsia toteutuksia. LlsSksi vastaanottimen tol- 
minnan ohjaamiseen voidaan kayttaa ohjausellnta, edullisesti mlkroprb- 
sessoria tai vastaavaa. 

10 Oheisessa kuvassa 7 on esitetty vieia keksinnon eraan edullisen 
suoritusmuodon mukainen elektronilkkalaite 24, joka kasitlaa 
langattoman vlestimen ja sljainnlnmaaritysvastaanottimeh. 
Ensimmaista antennia 13 kSytetaan paikannussatelliittien lahettamSn 
signaalin vastaanottami&edn. Vastaanotettu signaali Johdetaan 

15 enslmmalseen radlo-osaan 14, jossa signaali muunnetaan 
valltaaj'uudelle ja digitoidaan. Ensimmalnen radio-osa kSslttaa mm. 
kuvan 1 vastaanottimen muunninlohkon 2, digitointilohkon 3 ja 
kertojalohkon 4. Digltoitu signaali. joka tassa valheessa kasittaa 
edullisesti I- ja Q-komponentit, johdetaan digitaaliseen signaalikasittely- 

20 yksikkoGn 15, jossa suoritetaan mm. naytevektoreiden muodostus 
naytevektoreiden muodostuslohkossa 12. Naytteita tallennetaan 
ensimmalsiin muistivalineislln 16, jotka kasittavat mm. 
luku/kirjoitusmuistia seka edullisesti myos lukumuiatia ja/tai 
haihtumatonta luku/kirjoitusmulstia digitaalisen signaalinkasittely- 

25 yksikon 15 ohjelmakoodin tallennusta varten. Tassa suoritusmuodossa 
on signaalinkasittely-yksikossa 15 toteutettu mm. tahdistuslohkon 6 
toimintoja, kuten korrelaatiofunktiomatriisin C^^jt muodostus esim. 

Fourier-muuntimilla FF=T1. FFT2, .... FFTN ja/tai sovitetullla 
suodattimilia. Digitaalisessa signaalinkasittely-yksikossa 15 
30 muodostetaan edullisesti myds koherentt! etsintamatriisi Ax.k 1^ 

suoritetaan epakoherentti summausvaihe. Digitaalinen 
signaalinkasittely-yksikko 15 valittaa tiedon lasketusta vaihe-erosta ja 
taajuuspoikkeamasta suoritlnlohkoon 17, joka kasittaa esim. 
mikroprosessorin ja l/O-logiikkaa. Suoritinlohko 17 ohjaa 
35 skannauslohkoa 8 seka ensimmaista kytkinta 9. Seurantalohko 1 1 on 
edullisesti ainakin osittain toteutettu suoritiniohkon ohjelmakaskyina. 
Suoritinlohkon 17 datamuistina ja ohjelmamuistina kaytetaan toisia 
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muistivalineita 18, On selvaa, etta ensimmaiset muistivalineet 16 ja 
toiset muistivalineet 18 voivat kasittaa myos yhteista muistia. 
Sijainninmaaritysinformaatlota voidaan esittaa kayttajalle naytOlla 19, 

5 Suoritmlohkon 17 sovellusohjelmistossa on toteutettu my6s 
langattoman vlestimen toirmintoja. Nayttoa 19 voidaan talloin kayttaa 
sinansa tunnetusti myos mm. puheluinformaation esittamiseen. 
Nappaimistolla 20 kayttaja voi ohjata sijainnlmaarltysvastaanotinta seka 
langatonta viestinta. Koodekin 21 avulla suorit©«aan audiosignaaiien 

10 koodaaminen ja dekoodaaminen. Kuvassa 7 on esitetty viela ; 
langattoman viestimen radio-osa 22 la tolipi^n anteoni 2^ 

Nyt esilla olevaa keksintda ei ole rajoitettu ainoastaan edella esitGttyihin 
suoritusmuotolhin, vaan sita voidaan muunnella oheisten patenttlvaati- 
15 musten puitteissa. 
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Patenttivaatimukset: 

1. Menetelma vastaanottimen (1) tahdistamiseksi lahetettyyn koodimo- 
duloltuun hajaspektrisignaaliin, jossa menetelmassa kaytetaan ainakin 
5 yhta vertallukoodia (r(x)). joka vastaa jotakin moduloinnissa kaytettya 
koodia, selvitetaan lahetetyn signaalin taajuussiirtyma ja moduloinnissa 
kaytetyn koodin koodivaihe, lunnettu siita, etta menetelmassa 
suoritetaan ainakin seuraavia vaiheita: 

- naytevektoreiden muodostusvaihe, jossa vastaanotetusta slgnaallsta 
10 otetaan nSytteita (101) naytevektoreiden (Pk(1), Pk(2)...Pk(N)) 

m uodostamise ksi . 

- korrelaatiovaihe (102, 103. 104. 105) 
korrelaatlofunktiomatrlisin (Cx,k) muodostamiseksi naytevektoreiden 

(Pk(1), p,<(2).-Pk(N)) perusteella, 
15 - analyysivalhe (106, 107) koherentin etsintamatriisin (Ax,k) 

muodostamiseksi korrelaatlofunktiomatrlisin (Cx.k) transpoosin 
perusteella, ja 

- selvitysvailie (110), jossa selvitetaan taajuussiirtyma ja koodivaihe 
kayttamalla malnitun koherentin etsintamatriisin (Ax,k) a'ki- 

20 oiden ( a^^^f, (i , j) ) an/oja. 

2. Patenttivaatimuksen 1 mukainen menetelma, tunnettu slita, etta 
menetelmassa maaritetaan tutklttava taajuusalue, jaetaan maaritetty 
taajuusalue kahteen tai useampaan osaan, jolloin kullekin osalle 

25 suoritetaan mainltut naytevektoreiden muodostusvaihe, 
korrelaatiovaihe, ja analyysivalhe, tallennetaan koherentin 
etsintamatriisin (Ax,k) alkioiden (a^ jt(i.7)) an^oja. jolloin selvitysvalhe 

suoritetaan maarltetyn taajuusalueen tutkimisen jalkeen. ja 
taajuussiirtyma ja koodivaihe maaritetaan suurimman arvon 
30 perusteella. 

3. Patenttivaatimuksen 1 tai 2 mukainen menetelma, tunnettu siita, 
etta menetelmassa maaritetaan kynnysarvo, jolloin selvltysvaiheessa 
kaytetaan taajuussiirtyman ja koodlvalheen selvittamiseksi malnitun 

35 kynnysarvon ylittavia koherentin etsintamatriisin (Ax.k) 
alkloiden (o;, tO'.j)) suumusan/oja. 
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4. Patenttivaatimuksen 1 , 2 tai 3 mukainen menetelma, tunnettu silta, 
etta menetelmassa toistetaan (108) mainittuja nSytevektoreiden 
muodostusvaihetta, korrelaatiovaihetta ja analyysivaihetta ainakln 
kahden koherentin etsintamatriisin (Ax,k) muodostamiseksi, ja etta 

5 menetelmassa lisaksi suoritetaan summausvaihe (109). jossa 
muodostetaan epakoherentti etsintamatriisi (S^) summaamalla 
epakoherentisti kullakin toistokerralla muodostetut koherentin 
etsintamatriisin (Ax.k) vastlnalkloiden arvot (ax,/f ('.))). mainltussa 
selvitysvailieessa kaytetaan mainittua epakolierenttia etsintamatrlisla 

10 (Sji) taajuussilrtyman ja koodivaiheen maarityksessa. 

5. Patenttivaatimuksen 4 mukainen menetelma. tunnettu siita, etta 
menetelmassa maaritetasn tutkittava taajuusalue, jaetaan maaritetty 
taajuusalue kahteen tai useampaan osaan, jolloin kullekin osalle 

15 suoritetaan main'rtut naytevektoreiden muodostusvaihe, 
korrelaatiovaihe, analyysivaihe, ja summausvailie. seka tallennetaan 
epakolierentin etsintamatriisin {s^) alkioiden an/oja, jolloin selvitysvaiine 

suoritetaan maaritetyn taajuusalueen tutkimlsen jaikeen, ja 
taajuussilrtyma ja koodivailie maaritetaan suurimman an^on 
20 pemsteella. 

6. Patenttivaatimuksen 4 tai 5 mukainen menetelma, tunnettu siita, 
etta menetelmassa maaritetaan kynnysarvo, jolloin selvitysvaiheessa 
kaytetaan taajuussiirtyman ja koodivailieen selvittamlseksi mainitun 

25 kynnysarvon ylittavia epakoherentin etsintamatriisin (5;^) alkioiden 

suunjusarvoja. 

7. Jorikin patenttivaatimuksen 1 — 6 mukainen menetelma, tunnettu 
slita. etta kon'elaatiofunktlomatriisi (C^^k) muodostetaan sovltetuilla 

30 suodattimilla (26, 27, 28). 

8. Jonkin patenttivaatimuksen 1—7 mukainen menetelma, tunnettu 
siita, etta 

- korrelaatiovaiheessa muodostetaan mainittua vertailukoodia (r(x)) 
35 vastaavan takaperoisen koodin aika-taajuusmuunnos (R(x)), ja 

kunkin naytevektorin (Pk(l), Pk(2)...pk(N)) aika- 

taajuusmuunnos (Pk(i)), suoritetaan mainitun takaperoisen koodin 
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aika-taajuusmuunnoksen (R(x)) ja kunkin naytevektb- 
rin (Pk(l), Pk(2). .Pk(N)) aika-taajuusmuunnoksen (PrO)) valineh 
kertolasku, suoritetaan kultekin kertolaskutulokselie (Mx,k(i)) taajuus- 
aikamuunnos, ja muodostetaan kertolaskutulosten (Mx.k(l)) taajuus- 
aikamuunnoksista (m^k(0) korrelaatlofunktlomatrlisi ( Cx.k ) Ja etta 
- analyysivaiheessa koherentti etsintamatriisi (Ax.k) muodostetaan 

suorittamalla aika-taajuusmuunnos mainitun korrelaatlofunktiomat- 
riisin (Cx.k) tfaispoosin riveille. 



10 9. Jonkin patenttivaatimuksen 1 — 7 mukainen menetelma, tunnettu 
siit§, etta 

- korrelaatiovaiheessa muodostetaan mainitun vertallukoodin (r(x)) 
aika-taajuusmuunnos (R(x)), ja kutakin nSytevek- 
toria(pK(1), pk(2)...Pk(N)) vastaavan takaperoisen koodin aika- 

15 taajuusmuunnos suoritetaan mainitun vertallukoodin (r(x)) 

aika-taajuusmuunnoksen (R(x)) ja kunkin n§ytevekto- 
rin (pk(1), Pk(2)...Pk(N))_ takaperc^sen koodin aika- 
taajuusmuunnoksen (f^(/)) vSlinen kertolasku, suoritetaan kullekin 

kertolaskutulokselie (Mx.k(i)) taajuus-aikdmuunnos, ja muodostetaan 
20 kertolaskutulosten (Mx,k(i)) taajuus-alkamuunnoksista (mx.k(i)) 

korrelaatiofunktiomatriisi (Cx,k ). Ja eWS 

- analyysivaiheessa koherentti etsintamatriisi (Ax.k) muodostetaan 

suorittamalla aika-taajuusmuunnos mainitun korrelaatiofunktiomat- 
rilsin (Cx,k) transpoosin riveille. 

25 

to. Jonkin patenttivaatimuksen 1 — 1 mukainen menetelma, tunnettu 
slits, etta 

- korrelaatiovaiheessa muodostetaan mainitun vertallukoodin (r(x)) 
aika-taajuusmuunnos (R(x)), ja kunkin n&ytevek- 

30 torin (Pk(1). Pk(2)...Pk(N)) aika-taajuusmuunnos (Pi<(/)), muodostetaan 

vertailukoodin (r(x)) aika-taajuusmuunnoksesta (R(x)) 

kompleksikonjugaatti {R*(x)), suoritetaan mainitun 

kompleksikonjugaatin ja kunkin ndytevektorin (Pk(1). Pk(2)...pk(N)) 
aika-taajuusmuunnoksen (FV(/)) valinen kertolasku, suoritetaan 

35 kullekin kertolaskutulokselie (Mx,k(i)) taajuus-aikamuunnos. ja 
muodostetaan kertolaskutulosten (Mx,k(i)) taajuus- 

aikamuunnoksista (mx.k<i)) korrelaatiofunktiomatriisi (Cx,k).ja etta 
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- analyysivaiheessa koherentti Gtsintamatriisi (A^^k) rnuodostetaan 

suorittamalla aika-taajuusmuunnos mainitun korrelaatiofunktiomat- 
riisin (Cx,k) transpoosin riveille. 

5 11. Jonkin patenttivaatimuksen 1—7 mukainen menetelma, tunnettu 
siita, etta 

- korrelaatiovaiheessa rnuodostetaan mainitun vGrtailukoodin (r(x)) 
aika-taajuusmuunnos (R(x)), ja kunkin naytevek- 
torin (pk(1), Pk(2)..,Pk(N)) aika-taajuusmuunnos (Rc(/)), rnuodostetaan 

1 0 kustakin naytevektorin (Pk(1 ), Pk(2)...pk(N)) aika- 

taajuusmuunnoksesta (F\(/)) kompleksikonjugaatti (P*k(i) ). 

suoritetaan kunkin mainitun kompleksikbnjugaatin ja 
verlailukoodin (r(x)) aika-taajuusmuunnoksen (R(x)) vallnen 
kertolasku, suoritetaan kullekin kertolaskutulokselle (iVlxk(i)) taajuus- 
15 aikamuunnos, ja rnuodostetaan kertolaskutulosten (Mx,k(i)) taajuus- 
aikamuunnoksista (mx,k(i)) korrelaatiofunktiomatriisi (Cx,k)JS 

-analyysivaiheessa koherentti etsintamatriisi (Ax,k) rnuodostetaan 

suorittamalla aika-taajuusmuunnos mainitun korrelaatiofunktiomat- 
riisin (Cx,k) transpoosin riveille, 

20 

12. Jonkin patenttivaatimuksen 8 — 11 mukainen menetelma, tunnettu 
sitta, etta malnlttuna aika-taajtlusmuunnoksena kdytetaan Fourler- 
muunnosta, jolloin mainrttuna taajuus-aikamuunnoksena kaytetaah 
kaanteista Fourier-muunnosta. 

25 

13. Jonkin patenttivaatimuksen 1 — 12 mukainen menetelma, tunnettu 
siita, etta malnitussa selvitysvaiheessa (110) kaytetaan koherentin 
etsintamatriisin (Ax,k) alkioiden (a^jt(w)) suuruusarvoja. 

30 14. Patenttivaatimuksen 13 mukainen menetelma, tunnettu siita, etta 

mainitussa selvitysvaiheessa (1 10) kaytetaan koherentin 
etsintamatriisin (A^y^) alkioiden {a^ i,{ij)) suuruusan/ojen toiseen 

potenssiin korotetluja arvoja. 

35 15. Jonkin patenttivaatimuksen 1 — 14 mukainen menetelma, jossa 
selvitetaan vastaanottimen (1) sijarnti ja tallennetaan sljaintitteto. 
tunnettu silta, etta menetelmassa maadtetaan tutkittava taajuusaiue, ja 
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valitaan mainitulta taajuusalueelta taajuussiirtyman tutkimisen 
aloittamiseksi aloitustaajuus kayttamalla vastaanottimeen (1) 
tallennettua sijaintitietoa. 

5 16. Jonkin patenttivaatimuksen 1— 15 mukainen menetelma. jossa 
vastaanottimen (1) sijainti ja sijaintitieto selvitetaan kahden tai 
useamman satgllUtin lahettamast^ hajaspektrimoduloidusta signaalista, 
ja jossa menetelmassa kullekin mainitulle satelliitille valitaan yksiloiva 
koodi, tunnettu silta, etta menetelmassa suoritetaan mainittuja 

1 0 naytevektoreiden muodostusvaihetta, korrelaatiovaihetta, 

analyysivaihetta, ja selvitysvaihetta kunkin mainitun satelliitin 
lahettamalle signaalille. 

17. Vastaanotin (1), joka kasittaa ainakin tahdistusvalineet 
15 vastaanottimen (1) tahdistamiseksi lahetettyyn koodimoduloituun 
hajaspektrisignaaliin, jossa vastaanottimessa (1) on valineet (16) 
ainakin yhden vertailukoodin (r(x)) kayttamiseksi tahdistamisen 
yhteydessa, joka vertailukoodi (r(x)) vastaa jotakin modalolnnissa kay- 
tettya koodia, ja valineet (15) lahetetyn signaalln taajuussiirtyman ja 
20 moduloinnissa kaytetyn koodin koodivaiheen selvittamiseksi, tunnettu 
siita, etta vastaanotin (1) kasittaa lisaksi: 

- naytevektoreiden muodostusvallneet (12. 16) 
naytevektoreiden (pk(1), pk(2).»pK(N)) nnuodostamiseksi 
vastaanotetusta signaalista, 

25 - korrelaatiovaiineet (6. FFT1 , FFT2-..FFTN) 

korrelaatiofunktiomatriisin (Cx,k) muodostamlseksi naytevektoreiden 

(Pk(1), Pk(2)...pk(N)) perusteella, 

- analyysivalineet (6, 12) kolierentin etsintamatrlisin (Ax.k) 

nnuodostamiseksi korrelaatiofunktiomatriisin (Cx,k) perusteella, ja 
30 - selvitysvalineet (7, 17) taajuussiirtyman ja koodivaiheen 
selvittamiseksi kSyttamaiia mainittuja koherentin etsintamatriisin 

(Ax.k) (^jc,fc(^j)) an/oja. 

18, Patenttivaatimuksen 17 mukainen vastaanotin (1), tunnettu siita> 
35 etta se kasittaa valineet (17) tutkittavan taajuusalueen maarittamlseksc, 
valineet (5, 8, 9) maaritetyn taajuusalueen jakamiseksi kahteen tai 
useampaan osaan, jolloin kullekin osalle on jarjestetty suoritettavaksi 
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mainitut naytevektoreiden muodostus, korrelaatiofunktlomatrlisin (Cx,j^) 
muodostus, ja koherentin etslntamatnlsin (Ax,k) muodostus, vaiineet 
(16) koherentin etsintamatriisin (Ax,k) alkioiden (aj^jt(^y)) arvojen 
tallentamiseksi, ja valineet (15) taajuussiirtyman ja koodlvalheen 
5 maarlttamiseksi suurimman arvon perusteella. 

19. Patenttivaatimuksen 17 tal 18 mukainen vastaanotin (1). tunnettu 
siita, etta se kasittaa valineet (15) kynnysarvon maarittamiseksi, ja 
valineet (15) malnltun kynnysarvon ja mainittujen koherentin 

10 etsintamatriisin (Ax.k) arvojen vertaamiseksi taajuussiirtyman ja 

koodivaiheen selvittamiseksi. 

20. Patenttivaatimuksen 17, 18 tal 19 mukainen vastaanotin (1), 
tunnettu siita, etta malnittuja naytevektoreiden muodostusta, 

15 korrelaatiofunktiomatriisin (Cx.k) muodostusta Ja kolierentin 
etsintamatriisin (Ax,k) muodostusta on jarjestetty toistettavaksl 
koiierentin etsintamatriisin (Ax,k) muodostamlseksi ainakin kaksl 

kertaa, ja etta vastaanotin (1) kasittaa lisaksl summausvalineet 
epakoherentin etsintamatriisin {s^) muodostamlseksi summaamalla 
20 kullakin tolstokerraila muodostetut koherentin etsintamatriisin (Ax^k) 
vastinalkiolden arvot (a^fj^Cw))! jolloln mainitussa selvrtysvaiheessa 
kaytetaan mainlttua epakoherenttla etsintamatriisia {s^) 
taajuussiirtyman ja koodivaiheen maarityksessa. 

25 21- Patenttivaatimuksen 17 mukainen vastaanotin (1), tunnettu siita, 
etta se kasittaa valineet (17) tutkittavan taajuusalueen maarlttdmiseksi, 
valineet (5, 8, 9) maaritetyn taajuusalueen jakamiseksi kahteen tai 
useampaan osaan, jolloin kullekin osalle on jarjestetty suoritettavaksi 
mainitut naytevektoreiden muodostus, korrelaatiofunktiomatriisin (Cx^k ) 

30 muodostus. koherentin etsintamatriisin ( Ax,k ) muodostus, ja koherentin 
etsintamatriisin (Ax,k) summaus epakoherenttiln etsintamatrilsiin (5;^), 
valineet (16) epakoherentin etsintamatriisin (s^) aikioiden arvojen 
tallentamiseksi, Ja valineet (15) taajuussiirtyman ja koodivaiheen 
maarittamiseksi suurimman arvon perusteella. 

35 

22. Patenttivaatimuksen 21 mukainen vastaanotin (1), tunnettu siita, 
etta se kasittaa valineet (15) kynnysarvon maarittamiseksi, ja valineet 
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(15) mainitun kynnysarvon ja malnitun Gpakoherentin etsintamatriisin 
{S^) arvojen vertaamiseksi taajuussiirtyman ja koodivaiheen 

selvittamiseksi. 



5 23. Jonkin patenttivaatimuksen 17 — 22 mukainen vastaanot'm (1), 
tunnettu siita, etta se kasittaa sovitettuja suodattimia (26, 27, 28) 
korrelaatiof unktiomatriisin ( Cx^^ ) muodostamiseksi. 

24. Jonkin patenttivaatimuksen 17 — ^23 mukainen vastaanotin (1), 
10 tunnettu siita, etta: 

- korrelaatiovalineet kasitlavat val!neet(l6) mainittua 
vertailukoodia (r(x)) vastaavan takaperolsen koodin aika- 
taajuusmuunnoksen (R(x)) ja kunkin nayteyek- 

\ov\w (pk(1 ), pk(2)...Pk(N)) aika-taajuusmuunnoksen (Pk{i)) ' i| 

15 muodostamiseksi, v^lmeet (15) mainitun takaperoisen koodin aika-^^^^^^^^^^^^^^ 0^^^ 
taajuusmuunnoksen ( R(x)) ja kunkin naytevekto- 

rin (pk(1), pk(2)...Pk(N)) aika-taajuusmuunnoksen (Pk(i)) vSIisen 
kertolaskun suorittamiseksi, valineet (15) kullekin 
kertolaskutulokselle (Mx,k(i)) taajuus-aikamuunnoksen 
20 suorittamiseksi, ja valineet (15) korrelaatiofunktiomatriisin (C^^k) 

muodostamiseksi kertolaskutulosten (Mx.k(i)) taajuus- 

alkamuunnoksista (mx,k(i)). j^ etta 

- anaiyysivalineet kasittavat valineet (15) koherentin etsintamatriisin 
(Ax,k) muodostamiseksi suorrttamalla aika-taajuusmuunnos mainitun 

25 korrelaatiofunktiomatriisin (Cx,k) transpoosin riveille. 

25. Jonkin patenttivaatimuksen 17 — ^23 mukainen vastaanotin (1), 
tunnettu sllta, etta: 

- korrelaatiovalineet kasittavat valineet (16) mainitun ^ 
30 vertailukoodin (r(x)) aika-taajuusmuunnoksen (R(x)) ja kutakin 

uB)f\ey&V}iox\^ (pk(1)i pk(2)...pk(N)) vastaavan takaperolsen koodin 
aika-taajuusmuunnoksen (^(/)) muodostamiseksi, valineet (15) 

mainitun vertailukoodin (r(x)) aika-taajuusmuunnoksen (R(x)) ja 
kunkin naytevektorin (pk(1), pk(2)-..pk(N)) takaperolsen koodin aika- 
35 taajuusmuunnoksen (Hc(')) valisen kertolaskun suorittamiseksi, 

valineet (15) kullekin kertolaskutulokselle (M;(,k(i)) taajuus- 
aikamuunnoksen suorittamiseksi » ja vSilineet (15) 
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korrelaatiofunktiomatriisin (Cx.k) muodostamiseksl 
kertolaskutulosten (Mx.k(0) taajuus-aikamuunnoksista (mx.k(i)). ja eM 

- analyysivalineet kasittavat vSlineet (15) koherentln etsintamatriisin 
(Axk) muodostamiseksi suorittamalla aika-taajuusmuunnos malnitun 

5 korrelaatiofunktlomatriisin (Cx,k) transpoosin riveille/ 

26. Jonkin patsnttivaatlmuksen 17 — 23 mukainen vastaanotin (1). 
tunnettu siita, etta: 

- korrelaatlovalineet kasittavat vaiineet(16) mainitun 
10 vertailukoodln (r(x)) aika-taajuusmuunnoksen (R(x)) 

kompleksikonjugaatin (R*(x)) ja kunkin naytevek- 

torin (Pk(l ), pk(2)...Pk(N)) alka-taajuusmuunnoksen ( ) 

muodostamiseksi, vallneet (15) mainitun vertailukoodln (r(x)) aika- 
taajuusmuunnoksen kompleksikonjugaatin (R*(x)) |a kunkin 
1 5 naytevektorin (Pk(1 ). Pk(2)..-Pk(N)) aika-taajuusmuunnoksen ( ) 

valisen kertolaskun suorittamiseksi, vallneet (15) kullekin 
kertolaskutulokselle (Mx.k(i)) taajuus-alkamuunnoksen 
suorittamiseksi, ja vallneet (15) korrelaatiofunktiomatriisin (Cx,k) 

muodostamiseksi kertolaskutulosten (Mx.k(i)) taajuus- 

20 alkamuunnoksista (mx,k(i)), ja etta 

- analyysivalineet kasittavat vallneet (15) koherentln etsintamatriisin 
(Ax,k) muodostamiseksi suorittamalla aika-taajuusmuunnos mainitun 

korrelaatiofunktlomatriisin (C^.k) transpoosin riveille. 

25 27. Jonkin patenttivaatimuksen 17 — 23 mukainen vastaanotin (1), 
tunnettu slita, ettS: 

- korrelaatlovalineet kasittavat vallneet (16) kunkin naytevek- 
torin (Pk(i), Pk(2)...pk(N)) aika-taajuusmuunnoksen (Fk(0) 

kompleksikonjugaatin (P*k(i)) ja mainitun vertallukoodin (f(x)) aika- 

30 taajuusmuunnoksen (R(x)) muodostamiseksi, vallneet (15) mainitun 
vertailukoodln (r(x)) aika-taajuusrhuunnoksen (R(x)) ja kunkin 
naytevektorin (Pk(1). Pk(2)...pk(N)) alka- 
taajuusmuunnoksen kompleksikonjugaatin (P*k(i)) valisen 

kertolaskun suorittamiseksi, vSiineet (15) kullekin 
35 kertolaskutulokselle (Mx.k(i)) taajuus-aikamuunnoksen 
suorittamiseksi, ja vallneet (15) korrelaatiofunktlomatriisin (Cx,k) 
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muodostamiseksi kertolaskutulosten (Mx.k(i)) taajuus- 

aikamuunnoksista (nr^,K(i)) etta 
- analyysivallneet kasittavat valineet (15) koherentin etsintamatriisin 
(Axk) muodostamiseksi suorittamalla aika-taajuusmuunnos mainitun 
5 korrelaatiofunktiomatrilsin (Cx.k ) transpoosin riveille, 

28. Jonkin patenttivaatimuksen 24 — 27 mukainen vastaanotin (1), 
tunnettu siita, etta mainittuna aika-taajuusmuunnoksena on kaytetty 
Fourier-muunnosta, jolloin mainittuna taajuus-aikamuunnoksena on 

1 0 kaytetty kaanteista Fourier-muunnosta. 

29. Jonkin patenttivaatimuksen 17— 28 mukainen vastaanotin (1), 
tunnettu siita. etta selvitysvalineissa on kaytetty koherentin 
etsintamatriisin (Ax,k) alkioiden {a^^ti},j)) suuruusarvoja, kuten 

1 5 itseisarvoja. 

30. Jonkin patenttivaatimuksen 17—29 mukainen vastaanotin (1). 
tunnettu siita, etta selvitysvalineissa on kaytetty kolierentin 
etsintamatriisin (Ax,k) alkioiden {a^^f,{ij)) suuruusarvojen toiseen 

20 potenssiin korotettuja arvoja. 

31. Jonkin patenttivaatimuksen 17 — 30 mukainen vastaanotin (1), joka 
kasittaa valineet (14, 15, 17) vastaanottimen (1) sljainnin selvittAmiseksi 
ja valineet (18) sijaintitiedon tallentamiseksi, tunnettu siit&, etta 

25 vastaanotin (1) kasittaa lisaksi valineet (17) tutkittavan taajuusalueen 
maarittamiseksi, ja valineet (17, 22, 23) aloitustaajuuden valitsemiseksl 
main'rtulta taajuusalueelta vastaanottimeen (1) tallennetun sijaintitiedon 
perusteella. 

30 32. Elektroniikkalaite (24), joka kasittaa ainakin 
sijainninmaaritysvastaanottimen (1), tahdistusvalineet 
sijainninmaaritysvastaanottimen (1) tahdistamiseksi lahetettyyn 
koodimoduloituun hajaspektrisignaaliin, jossa 

sijainninmaaritysvastaanottimessa (1) on valineet (16) ainakin yhden 

35 vertallukoodin (r(x)) kdyttSmiseksi tahdistamisen yhteydessa, joka 
vertailukoodi (r(x)) vastaa jotakin moduloinnissa kaytettya koodia, ja 
joka elektroniikkalaite kasittaa valineet (15) lahetetyn signaalin 
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taajuusslirtyman ja moduloinnissa kaytetyn koodin koodivalheen 
selvlttamiseksi. tunnettu siita. etta elektroniikkalalte (24) kasittaa 
lisaksi: 

- naytevektoreiden muodostusvallneet (12, 16) 
5 naytevektoreiden (pk(1),pk(2)...Pk(N)) muodostamiseksi 

vastaanotetusta signaalista, 

- korrelaatiovalineet (6, FFT1, FFT2...FFTN) 
korrelaatiofunktiomatriisin (Cx,k) muodostamiseksi naytevektoreiden 

(Pk(1), Pk(2)...Pk(N)) perusteella, 
10 - analyysivaiineet(6, 12) koherentin etsintamatriisin (Ax.k) 
muodostamiseksi konrelaatiofunktiomatriisin (Cx,k) perusteella, ja 

- selvitysvallneet(7, 17) taajuussiirtymSn ja koodivaiiieen 
selvittamiseksl kayttamaila mainittuja koherenttien etsintamatriisin 
( Ax.k ) alkioiden {ax,k{iJ)) arvoja. 

15 

33. Patenttivaatimuksen 32 mukalnen elektroniikkalaite (24), joka 
kasittaa valineet(l4, 15, 17) elektronllkkalaWeen (24) sijalnnin 
seivlttamiseksi ja vaiineet (18) sijaintitiedon tallentamiseksl. tunnettu 
siita, etta elektroniikkalaite (24) kasittaa llsaksl vaiineet (17) tutkittavah 

20 taajuusalueen maarlttamiseksl. ja vaiineet (17, 22, 23) aloitustaajuuden 
valitsemiseksi mainitulta taajuusalueetta vastaanottimeen (1) 
tallennetun sijaintitiedon perusteella. 

34. Patenttivaatimuksen 32 tai 33 mukalnen elektronrikkalaite (24), 
25 tunnettu siita, etta se kasittaa lisaksi vaiineet matkaviestintoimintojen 

suorittamiseksi- 

35. Patenttivaatimuksen 34 mukainen elektroniikkalaite (24), tunnettu 
siita, etta se kasittaa vaiineet (22, 23) tiedonsiirtoytiteyden 

30 muodostamiseksi matkaviestinverkkoon. jolloin vaiineet 
elektroniikkalaitteen (24) sijainnin selvittamiseksi kaslttavat 
vaiineet (22, 23) sijainninmaarityksessa kaytettavan tiedon hakemiseksi 
matkaviestinverkosta. 
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(57) Tiivistelma : 

Keksinnon kohteena on menetelma vastaanottimen (1) 
tahdistamiseksi lahetettyyn koodimoduloituun hajaspekt- 
risignaaliin. Menetelmassa kaytetaan ainakin yhta vertai- 
lukoodla (r(x)), joka vastaa jotakin moduloinnissa kaytet- 
tya koodia. IVIenetelm&ssa selvitetaan lahetetyn signaa- 
lin taajuussiirtyma \a moduloinnissa kaytetyn koodin 
koodivaihe. MenetelmSssa suoritetaan lisaksi ainakin 
seuraavia vailieita: 

- naytevektoreiden muodostusvaihe, jossa vastaanote- 
tusta signaalista ototaan naytteita (101) naytevektorei- 
den (pk(1). pk(2)...pk(N)) muodostamiseksi, 

- korrelaatiovaihe (102, 103, 104, 105) kon-elaatiofunk- 
tlomatriisin (Cx,k) muodostamiseksi naytevektoreiden 

(pk(1), pk(2)...pk(N)) perusteella, 

- anaiyysivailie (106, 107) kolierentin etsintamatriisin 
(Ax,k) muodostamiseksi kon-elaatiofunktiomatrilsin 
(Cx,k) franspoosin penjsteella, ja 

' selvitysvaihe (110). jossa selvitetaan taajuussiirtyma ja 
koodivaihe kayttamalla mainitun koherentin etsintamat- 
riisin ( Ax,K ) alkioiden (ax,*(».i)) arvoja. 



Fig. 1 
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Method for synchronizing a receiver, and a receiver 

The present invention relates to a method according to the preamble of 
Claim 1 for synchronizing a receiver, a receiver according to the preamble of 
Claim 1 7 and an electronic device according to the preamble of Claim 32. 

A well known positioning system is the GPS system (Global Positioning 
System), which comprises over 30 satellites, of which a maximum of 12 
satellites are simultaneously visible to a receiver. These satellites transmit 
various information, such as Ephemeris data and the satellite time by the 
clock. Normally, the receiver used in positioning determines its location by 
calculating the transit time to the receiver of a signal transmitted 
simultaneously from several satellites of the positioning system. For 
determining its location, the receiver must typically receive the signal of at 
least four visible satellites in order to be able to calculate the location. 

Each satellite of the GPS system in operation transmits a so-called L1 signal 
at the carrier frequency of 1575.42 MHz. This frequency is also denoted 
154fo, where fo=10.23 MHz. In addition, the satellites transmit an L2 signal at 
the carrier frequency of 1227.6 MHz, or 120fo. In the satellite, these signals 
are modulated with at least one pseudo-random sequence. Each satellite 
has a different pseudo-random sequence. After the modulation, the signal is 
a code modulated wideband signal. The modulation technique used makes it 
possible for the receiver to distinguish between signals transmitted by 
different satellites, although the carrier frequencies used in the transmission 
are essentially the same. This modulation technique is called Code Division 
Multiple Access (CDMA). The pseudo-random sequence used in each 
satellite for the modulation of the LI signal is, for instance, a so-called C/A 
code (Coarse/Acquisition code), an example of which is the Gold code. Each 
GPS satellite transmits a signal by using a unique C/A code. The codes are 
formed as the modulo 2 sum of two 1023-bit binary sequences. The first 
binary sequence G1 has been formed with the polynome X^°+X^+1, and the 
second binary sequence G2 has been formed by delaying the polynome 
X^°+X^+X®+X^+X^+X2+1 so that each satellite has a different delay. This 
arrangement enables generating different C/A codes with a similar code 
generator. The C/A codes are binary codes, the chipping rate of which in the 
GPS system is 1.023 MHz. The C/A code comprises 1023 chips, whereby 
the repeating time of the code is 1 ms. The carrier wave of the L1 signal is 



further modulated with navigation information at a bit rate of 50 bit/s. The 
navigation information comprises information about the "health", orbit and 
time of the satellite, for instance. 

The satellites monitor the condition of their equipment during operation. For 
example, the satellites can use so-called watch-dog functions to detect 
errors in the equipment and notify about them. Errors and malfunctions can 
be either momentary or long-term in nature. On the basis of the health 
information, some of the errors can be possibly compensated for, or the 
information transmitted by a malfunctioning satellite can be entirely ignored. 
In addition, in a situation in which the signal of more than four satellites can 
be received, the information received from different satellites can be 
emphasized in different ways on the basis of the health information. Thus it 
is possible to minimize errors that satellites which seem to be unreliable may 
cause in the measurements. 

In order to detect the signals of the satellites and to identify the satellites, the 
receiver must perform synchronization, whereby the receiver finds the signal 
of each satellite each time and tries to synchronize with this signal so that 
the data transmitted with the signal can be received and demodulated. 

In prior art receivers, the time needed for this synchronization depends on 
the strength of the received signal, among other things. Typically, the 
weaker the signal to be received is, the longer integration is needed in each 
element of the search space (correlation/frequency) in order to detect a 
possible signal. In prior art GPS receivers, which are typically designed for 
outdoor use, synchronization with the signals of the satellites takes some 
tens of seconds or a few minutes, if the received signal is relatively strong, in 
the range of -120 - -130 dBm. However, if the positioning is to take place 
indoors or In a place where the received signal is attenuated by buildings or 
other terrain obstacles, the synchronization time increases substantially. A 
typical signal strength indoors is approx. -150 dBm, whereby the time 
needed for the synchronization process with a sequential search would be 
hours. This time is clearly too long for performing a reasonable positioning. 
Figure 2a illustrates, by way of example, a CDMA signal transmitted from a 
satellite, and correspondingly Fig. 2b illustrates this transmitted signal when 
it reaches the receiver. The signal to be received has been attenuated 
substantially on the transmission link, and it contains accumulated noise. 



The positioning receiver must perform synchronization when, for instance, 
the receiver is switched on, and in a situation that the receiver has not been 
able to receive the signal of any satellite for a long time. A situation like this 
can arise easily in mobile devices, for example, because the device moves 
and the antenna of the device is not always in an optimal position In relation 
to the satellites, which weakens the strength of the signal coming to the 
receiver. In city areas, buildings also have an effect on the received signal, 
and in addition, the signal can proceed to the receiver along several different 
routes, such as straight from the satellite (line-of-sight) and as reflected from 
buildings. As a result of this multipath propagation, the same signal is 
received as several signals in different phases. 

The positioning arrangement has two main functions: 

1 . calculating the pseudo distance of the receiver to different GPS 
satellites, and 

2. determining the location of the receiver using the calculated 
pseudo distances and the location information of the satellites. 
The location information of the satellites at any given time can be 
calculated on the basis of the Ephemeris and time correction 
information received from the satellites. 

The distances to the satellites are called pseudo distances, because the 
time is not known accurately in the receiver. Then the determination of the 
location and time is repeated, until a sufficient accuracy of time and location 
has been achieved. Because the time is not known with an absolute 
accuracy, . the location and time must be determined by linearizing the 
equation group for each new iteration. 

The calculation of the pseudo distance can be carried out by measuring the 
mutual, apparent transit time delays of the signals of different satellites. After 
the receiver has become synchronized with the received signal, the 
information transmitted in the signal is clarified. 

Almost all known GPS receivers use correlation methods for calculating the 
distances. The pseudo random sequences of different satellites have been 
saved or are generated locally in the location determination receiver. A down 



conversion is performed on the received signal, wliereafter tlie receiver 
multiplies the received signal with the saved pseudo random sequence. The 
signal formed as the result of the multiplication is integrated or low-pass 
filtered, whereby the result tells if the received signal contained a signal 
transmitted by the satellite. The multiplication carried out in the receiver is 
repeated so that each time the phase of the pseudo random sequence 
saved in the receiver is shifted. The right phase is concluded from the 
correlation result preferably so that when the correlation result is the highest, 
the right phase has been found. Then the receiver is correctly synchronized 
with the received signal. 

After the synchronization with the code has been performed, the frequency 
is adjusted and phase locking is carried out. This correlation result also 
reveals the information transmitted in the GPS signal, which means that it is 
a demodulated signal. 

The above mentioned synchronization and frequency adjustment process 
must be repeated on each signal of a satellite which is received in the 
receiver. Therefore, this process consumes a lot of time especially in a 
situation that the received signals are weak. In order to make this process 
faster, several correlators are used in some prior art receivers, whereby 
more correlation peaks can be searched for simultaneously. Only by 
increasing the number of correlators it is not possible to make the 
synchronization and frequency adjustment process much faster in practical 
applications, because the number of correlators can not be increased to 
infinity. 

Some prior art GPS receivers use the FFT technique for determining the 
doppler shift of the received GPS signal in connection with ordinary 
correlators. In these receivers, correlation is used to reduce the bandwidth of 
the received signal to 10 kHz-30 kHz. This narrow-band signal is analyzed 
with FFT algorithms to determine the carrier frequency. 

A GPS receiver and a method for processing GPS signals have been 
presented in the international patent application WO 97/14057. In the main, 
the receiver presented in this specification comprises two separate 
receivers, of which the first receiver is intended for use in a situation that the 
strength of the received signal is sufficient, and the second receiver is 



intended for use in a situation that the strength of the received signal is not 
sufficient for a sufficiently accurate location determination when the first 
receiver is used. In this second receiver, the received signal is digitized and 
saved in memory media, whereby these saved signals are later processed in 
a digital signal processing unit. The digital signal processing unit performs 
convolution operations to the received, digitized signal. The purpose of 
these convolution operations is to calculate pseudo distances. The number 
of PM frames saved in the memory media is typically from one hundred to 
one thousand, which corresponds to a signal with the length of 100 ms to 
1 s. After this, a saved code corresponding to the code of the satellite to be 
examined is retrieved from the memory of the receiver to be used for 
analyzing the received signal. 

The doppler shift is also removed in the receiver. The amount of this doppler 
shift has been determined either with the first receiver or on the basis of 
information received from a base station of the GPS system. After this, a 
coherent summing of consecutive frames is performed. A fast Fourier 
transformation is performed on the data set received as the result of this 
summing. A multiplication is performed on the Fourier transformation result, 
using the complex conjugate of the Fourier transformation of the reference 
signal saved in the memory media as one multiplier. An inverse Fourier 
transformation is further performed on the result of this multiplication, 
whereby a number of correlation results are received. In this specification, 
correlation has been replaced by a Fourier transformation, whereby the 
number of calculations has been reduced. According to the specification, the 
method makes location determination 10 to 100 times faster as compared to 
the solutions known at the time of filing the specification. 

It is an objective of the present invention to provide a receiver, in which 
synchronization with the transmitted signal can be performed substantially 
faster than in the prior art receivers, even if the signals are weak. The 
invention is especially suitable for use in location determination receivers, 
but also in other receivers, preferably in CDMA receivers, in which the 
receiver must become synchronized with a spread spectrum signal. The 
invention is based on the idea that the frequency analysis properties of 
convolution and time-to-frequency transformation, such as the fast Fourier 
transformation, are utilized more efficiently, whereby it is possible to 
examine a greater amount of information simultaneously than in the prior art 



solutions, which use the Fourier transformation. The method according to 
the present invention is characterized in what Is set forth in the 
characterizing part of Claim 1. The receiver according to the present 
invention is characterized in what is set forth in the characterizing part of 
Claim 17. The electronic device according to the present invention is 
characterized in what is set forth In the characterizing part of Claim 32. 

The present invention provides considerable advantages as compared to the 
prior art methods and receivers. Especially In location determination 
receivers, the method according to the invention makes It possible for the 
receiver to become synchronized with the received signal faster. Then the 
location can be determined much faster even when the received signal is 
weak, which is the case Indoors or in city areas where there are obstacles 
that substantially attenuate the propagation path of the signal. The method 
according to the Invention provides a substantial Improvement of the 
synchronization speed as compared to the prior art methods. The receiver 
according to the invention can be realized with a relatively small number of 
components, and the total energy consumption can be kept reasonable, 
whereby the invention is especially suitable for use in portable devices. Thus 
the location determination receiver can also be Implemented in connection 
with a mobile station. 



In the following, the present invention will be described In more detail with 
reference to the accompanying drawings, in which 

Figure 1 Is a simplified block diagram of a receiver, in which the method 
according to the Invention can be applied. 

Figure 2a shows the signal transmitted from the satellite as a simplified 
diagram, 

Figure 2b shows an example of the simplified signal of Fig. 2a in a receiver. 

Figure 3 shows the correlation step of a method according to a preferred 
embodiment of the invention as a simplified diagram. 



Figure 4 shows the analysis step of a method according to a preferred 
embodiment of the invention as a simplified diagram. 



Figure 5 shows the summing step of a method according to a preferred 
embodiment of the invention as a simplified diagram, 

Figure 6 shows the inference step of a method according to a preferred 
embodiment of the invention as a simplified diagram, 

Figure 7 shows an electronic device according to a preferred embodiment 
of the invention as a simplified block diagram, and 

Figure 8 shows a prior art correlator, which can be used in connection 
with a method according to a preferred embodiment of the 
invention. 

In the receiver 1 shown in Fig 1 , the received signal is preferably converted 
to an intermediate frequency in the converter block 2. At this stage, as is 
well known, the signal comprises two components, I and Q, between which 
there is a phase difference of 90°. These analog signal components, which 
are converted to the intermediate frequency, are digitized in the digitizing 
block 3 and conducted to the multiplier block 4. In the multiplier block 4, the I 
and Q components of the digitized signal are multiplied with the signal 
generated by the Numerically Controlled Oscillator (NCO) 5. This signal of 
the numerically controlled oscillator 5 is intended to correct the frequency 
deviation caused by the doppler shift and the frequency error of the local 
oscillator (not shown) of the receiver 1. This numerically controlled oscillator 
5 is also used in connection with the method according to the invention, as 
will be described later in this specification. The signal formed by the 
multiplier block 4 is led to the synchronization block 6, in which the method 
according to the invention is applied. This synchronization block 6 tries to 
find the code phase and frequency deviation of the satellite to be used in 
functions after the synchronization. This will also be described later in this 
specification. The frequency of the numerically controlled oscillator 5 is 
regulated with the scanning block 8 when required, and the scanning block 8 
is controlled with the control block 7 during synchronization. The control 
block 7 controls the first switch 9 to switch either the signal formed by this 
scanning block 8 to the numerically controlled oscillator 5 during 
synchronization or the control signal formed by the code tracking block 1 1 to 
the numerically controlled oscillator 5 when synchronization has been 



achieved. The second switch 10 is used to control the operation of this 
tracking block 1 1 . This tracking block 1 1 forms a part of a code phase locked 
loop and a carrier phase locked loop, which are known as such (not shown). 

5 After switching the operating voltages on or in a situation that the receiver 1 
has not been able to receive a signal of the GPS satellites for a long time, a 
two-dimensioned search step is carried out in the receiver 1 for each 
satellite to which signal reception is targeted. In this two-dimensioned search 
step, the purpose is to determine the carrier frequency and code phase of 

10 each satellite. This carrier frequency is influenced by the doppler shift 
caused by the movement of the satellite and the inaccuracies of the local 
oscillator of the receiver. The frequency inaccuracy can be rather high, as 
much as ± 6 kHz, whereby the receiver must search from the frequency 
range of approx. 12 kHz in relation to the actual transmission frequency 

15 (LI = 1575.42 MHz). The receiver 1 does not know the exact code phase, 
either, whereby the receiver must also determine the code phase from 1023 
possible code phases. This results in a two-dimensioned search process, in 
which one dimension is a frequency deviation in the range of 12 kHz, and 
the other dimension is a code phase out of 1023 different code phases. In a 

20 method according to a preferred embodiment of the invention it is possible to 
examine a frequency range of approx. 500 Hz at a time, whereby the 
method is repeated 24 times when required, to cover the whole frequency 
range of 12 kHz to be examined. Obviously the values used in this 
specification only serve as examples which clarify the invention, but do not 

25 constitute a limitation thereof. The invention can also be applied in other 
than GPS systems, whereby the frequency values, code phases and the 
number of codes can vary. 

The operation of a method according to a preferred embodiment of the 
30 invention in a receiver 1 according to Figure 1 will be described in the 
following. In order to start synchronization, the scanning block 8 sets the 
frequency of the numerically controlled oscillator 5 so that the receiver 
advantageously receives the lowest frequencies of the frequency range, 
1575.414 MHz — 1575.4145 MHz. The receiver can also specify the starting 
35 frequency so that the receiver utilizes previously determined location 
information and/or almanac information, whereby the location determination 
can be made even faster. Samples of the received signal are saved in the 
sample vector formation block 12 preferably as complex sample vectors 



Pk(1), Pk(2)...Pk(N), each of which includes 1023 samples in this preferred 
embodiment. In this preferred embodiment, the sample saving frequency of 
the sample vector formation block 12 is essentially the same as the chipping 
rate of the chips, which is approx. 1 023 000 samples in a second. The 
sample vectors are continuous so that the next sample vector continues 
temporally after the previous sample vector, or the time difference between 
the last sample of the sample vector and the first sample of the next sample 
vector is essentially the same as the time difference between consecutive 
samples of the sample vector. These 1023 samples thus correspond to a 
signal of 1 ms, whereby in a time-to-frequency transformation the frequency 
range is 1 kHz, which can be partly utilized. The sample vector formation 
step is denoted with the reference number 101 in Fig. 3 in this specification. 

The number of sample vectors is preferably N, where N is advantageously a 
power of two. In addition, in a preferred embodiment of the invention the 
formation of the sample vectors Is repeated K times, as will be described 
later in this specification. In the following denotations, the subscript k is used 
for different times of repetition. When the value of the number N of sample 
vectors Pk(1), Pk(2)...pk(N) is determined in the GPS system, it must be taken 
into consideration that information has been modulated on to the signal at a 
bit rate of 50 bits/s as binary phase modulation. Another factor, which limits 
this number N of the sample vectors pk(1), Pk(2)...pk(N), is the frequency 
stability of the local oscillator of the receiver. 

In addition to the sample vector formation step, the synchronization method 
according to the invention also comprises a correlation step for forming a 
correlation function matrix. 

This correlation step can be partly performed already during sampling, or 
after N sample vectors Pk(1), Pk(2)...pk(N) have been formed. If the 
correlation step is performed in such a manner, for instance, that a time-to- 
frequency transformation, such as a Fast Fourier Transformation (FFT), is 
calculated for each sample vector after it has been saved, the same time-to- 
frequency transformer can be used in all N sample vectors 
Pk(1), Pk(2)...pk(N). But if the correlation step is performed after the saving of 
N sample vectors, a separate time-to-frequency transformer must be used 
for each sample vector, or the time-to-frequency transformations are 
performed for different sample vectors consecutively in the same time-to- 
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frequency transformer. Fig. 3 shows the correlation step of the method, in 
which a correlation function matrix C^.k -is formed of the sample vectors 

Pk(1), pk(2)...pk(N). Although the Fourier transformation will be mainly used 
as an example of a time-to-frequency transformation later in this 
5 specification, and an inverse Fourier transformation will be used as an 
example of an inverse transformation, or a frequency-to-time transformation, 
it is obvious that this invention is not limited to these examples only. 

A discrete Fourier transformation 102, most suitably a fast Fourier 
10 transformation, FFT, is performed on each sample vector Pk(1), Pk(2)...pk(N). 



This Is illustrated by the blocks FFT1, FFT2, FFTN in Fig. 3. The number 
15 of values used in calculations is preferably 1024, because then the discrete 
Fourier transformation can be implemented in practical applications much 
more efficiently (with an FFT algorithm) than when 1023 values are used. 
One way of doing this is to add an extra null as the 1024th element. This has 
a very small effect on the transformation result. 



The reference codes r(x) corresponding to the C/A code of all the satellites 
of the GPS system are preferably saved in the receiver, whereby x refers to 
the satellite identifier and is, for example, in the range 1 — 36. It is not 
necessary to save the reference codes, but they can also be generated in 
25 the receiver. The reference code, which is selected or generated in the 
correlation step, is the reference code of the satellite transmitting the signal 
with which the receiver is being synchronized at any given time. The 
reference code is temporally inversed. This inverse reference code, which is 
denoted with f(x) in Fig. 3, Is treated with a discrete Fourier transformation 

30 1 03, most suitably a fast Fourier transformation (FFT). 



The inverse reference code r(x) and/or its FFT transformation may have 

35 been saved in the memory media of the receiver in advance, or it is formed 
of the reference code r(x) in connection with synchronization. 




(1) 



20 



R(x) = FFT(f(x)) 



(2) 



A multiplication 104 between the Fourier transformation result Pk(i) of each 
sample vector Pk(i) and the Fourier transformation R(x) of the inverse 

reference code f{x) is performed next in the correlation step. 



Finally, an inverse Fourier transformation 105 is performed on the results of 
these multiplications, whereby the result Is a cross-correlation of the 
reference code r(x) and the received signal with all possible Integer delays 
(1023). 



This Is based on the fact that the Fourier transformation of the convolution of 
the time domain signals corresponds to the multiplication of Fourier- 
transformed signals, or time domain signals converted to the frequency 
domain. When an inverse reference code Is also used, it is possible to 
perform a fast discrete time correlation with the Fourier transformation. In 
this preferred example, the cross-correlation result then comprises 1023 
elements. The cross-correlation results mx,k(i) formed of these different 
sample vectors Pk(i) are used to form a correlation function matrix Cx,k . 

which the number of rows is the number N of sample vectors. 

It Is obvious that instead of temporal inversion of the reference code, it is 
possible to form inverse sample vectors pk(i) of the sample vectors Pk(l), 

whereby the reference code r(x) is used as direct and the sample vectors 
are used as inverse in the calculations presented above. In a preferred 
embodiment, it is not necessary to perform either of the inversions 
mentioned above, but the reference code r(x) and the sample vectors Pk(i) 
can be used as such. This is based on utilizing the properly of the 
correlation theorem, which expresses that the cross-correlation corr(Zi, Za) 
between two time discrete functions Zi, Z2 can be formed by means of a 
frequency-to-tlme conversion of functions converted to the frequency 
domain. This is expressed with the following formula: 



Mx.k(i) = R(x)-Pk(i) 



(3) 




corr(zi, Zz) IFFt(z/ • Z2) 
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In a preferred embodiment of the invention, it is tfius possible to form the 
cross-correlation between the sample vectors pk(i) and the reference code 

r(x) by performing a Fourier transformation (most suitably a fast Fourier 
5 transformation FFT) on the sample vectors p^^{\) and the reference code r(x) 

for forming a Fourier transformation Pk(i) of the sample vector Pk(i) and a 

Fourier transformation R(x) of the reference code r(x) , by forming a complex 

conjugate P*k(i) , of the Fourier transformation Pk(i) of the sample vector, by 

multiplying the formed complex conjugate P*k(i) with the Fourier 

10 transformation R(x) of the sample vector and by performing an Inverse 

Fourier transformation on the multiplication result. As an alternative, a 
complex conjugate R*(x) can be formed of the Fourier transfomriation R(x) 

of the reference code, whereby it is multiplied with the Fourier transformation 
Pk(i) of the sample vector, whereafter an inverse Fourier transformation is 

15 performed on the multiplication result. 

In this connection, it should be emphasized that the method presented in the 
previous paragraph for calculating the cross-correlation between the sample 
vectors Pk(i) and the reference code r(x) is derived from the basic properties 

20 of correlation and convolution and the close dependence between them, 
whereby the inversion of the function in the time domain is practically 
mathematically dependent on the formation of a complex conjugate in the 
frequency domain. This is dealt with in more detail e.g. In the specification: 
"Digital Signal Processing - A Practical Approach", Emmanuel C. Ifeachor 

25 and Barrie W. Jervis, Addison-Wesley Publishing Company 1993, ISBN 0- 
201 -5441 3-X, paragraph 4: "Correlation and Convolution", which is referred 
to in this connection. It should also be mentioned that as a whole, it is not 
important for the application of this invention which method is used for 
achieving the cross-correlation result. 

30 

The rows of a correlation function matrix Cx,k formed in the correlation step 

show the cross-correlation of the received signal and the reference code 
with different phase differences at the intervals of one millisecond. As a 
formula, the correlation function matrix can then be expressed as follows: 

35 
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mx.k(l)' 
mx.k(2) 



(6) 



mx.k(N) 



At the next or analysis phase, the transpose 1 06 of the correlation function 
matrix Cx,k is used, in which transpose the rows represent the signal 

samples in the time domain In the same manner as In a prior art correlator. 
Each row corresponds to a certain code phase difference between the 
received signal and the reference code. A Fourier transformation 107 is 
performed on each row of the transpose of this correlation function matrix 
Cx,k for forming a coherent search matrix Ax,k , whereby a frequency 

analysis can be performed for determining the real frequency shift. This is 
Illustrated by Fig. 4. 



In practical applications, it is not necessary to form a transposed matrix from 
a correlation function matrix In a separate step, but the elements of the 
saved correlation function matrix Cx.k are read from the memory 16 (Fig. 7) 

In a different direction, preferably by columns. 

The correlation function matrix Cx,k can also be formed e.g. by using 
matched filters in a manner known as such. For example, the patent 
specification US-5,373,531 discloses a correlator 25 equipped with matched 
filters (Fig. 8), which can be applied for forming a correlation function matrix 
Cx,k- The matched filter is formed of three shift registers so that the I- 

component of the received signal Is led to the first shift register 26, the Q- 
component of the received signal Is led to the second shift register 27, and a 
reference code corresponding to the code of the satellite Is led to the third 
shift register 28. The length of the shift registers 26, 27, 28 is preferably the 
same as the length of the sample vector, which is 1023 bits in the GPS 
system. When one sample vector is transferred to the shift register 26, 27, 
the correlation is performed as follows. An exclusive NOR operation is 
performed between the bits of the shift register 26 containing the I- 
component and the bits of the shift register 28 containing the reference code 




(7) 
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in the first reference block 29, and an exclusive NOR operation is performed 
between the bits of the shift register 27 containing the Q-component and the 
bits of the shift register 28 containing the reference code in the second 
reference block 30. The result of the operation is bit-specific correlation 
5 infornnation: if a bit of the received signal (examined by I and Q components) 
and a corresponding bit of the reference code match, the result is a binary 
value 1 , and if not, the result is a binary value 0. If the signal correlates 
completely with the reference code, all the bits have the value 1 . If none of 
the bits correlate, all the bits have the value 0. The values received as the 

10 result of the correlation are summed in the first summing block 31 (l-values) 
and in the second summing block 32 (Q-values), whereby the result is the 
correlation value mx,k(i) with a certain code delay in a certain time. After one 
comparison, the reference code is shifted by one, whereafter a new 
correlation result is formed. Each correlation result formed by the method 

15 described above corresponds to one element of a row of the correlation 
function matrix Cx,k (1023 in total). The steps mentioned above are 

repeated on each sample vector when required, whereby one row of the 
correlation function matrix Cx,k is obtained at each time of calculation. 

20 However, as was mentioned above, in a GPS system the signal is 
modulated with a signal of 50 bit/s, which limits the value of the number N in 
practical applications. Then the number N must be selected preferably so 
that the modulation does not have a substantial effect on the analysis. In 
addition, the optimum value of N is dependent on the selection of a window 

25 function for the Fourier transformation. If the numerical value is selected as, 
for instance, N = 32, the width of the noise band is In the range of 30 Hz, 
which is still a little too much for expressing signals with a strength of the 
order of -150 dBm in the receiver. Because of this, an optional incoherent 
summing step, in which the signal-to-noise ratio is improved, is performed in 

30 the synchronization block 6. 

In order to implement the incoherent summing step, the sample vector 
formation step, the correlation step and the analysis step described above 
are repeated 108 K times (Fig. 5). This number of repetitions K is preferably 
35 selected such that the signal-to-noise ratio can be improved sufficiently, but 
in a reasonable time. One coherent search matrix Ax^^ 's formed each time 

when the analysis step is performed, and an incoherent summing operation 
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is performed on it to form an incoherent search matrix . An incoherent 
search matrix is preferably formed as follows. To advantage, the quantity 
(absolute value) or other numerical value, such as the second power of the 
quantity value of the element, is calculated from the complex elements 
^xA''J) esich coherent search matrix Ax,k- The numerical values 
calculated from the equivalent elements of each incoherent search matrix 
are summed, or in other words, an addition of the matrixes is performed, 
which can be expressed with the formula: 



K 

k-1 



hx.k(l.1)| - |ax.k(1.N)| 
a^k(l023,1)| ... |a^,(l023.N)|_ 



(8) 



In practical applications, an incoherent search matrix can be formed in at 
least two ways. In the first alternative implementation, a coherent search 
matrix formed at each time of repetition is saved. After the required times of 
repetition, an incoherent search matrix is performed by summing the 
equivalent elements according to Formula 8. In this alternative 
implementation, memory is needed for saving the elements of all the 
coherent search matrixes. According to another alternative implementation, 
one coherent search matrix is calculated at first, and the values of this matrix 
are copied as the elements of the sum matrix. A coherent search matrix is 
formed at each following time of repetition, and the values of the matrix are 
summed with the equivalent elements of the incoherent search matrix. In 
this alternative, the summing of the equivalent elements is thus performed at 
each time of repetition. Only one coherent search matrix is thus saved, 
whereby less memory is needed than in the first alternative. 

After the required times of repetition have been performed, the values of the 
elements Sx(i,j) of this incoherent search matrix are examined in the 

determination step, and the objective is to find 110 a value which exceeds 
the threshold value set in advance and is clearly higher than the other 
values. If a value like this is found, it expresses the code phase difference 
and frequency deviation, because it is probably a signal transmitted by a 
satellite. If the signal is not transmitted by a satellite, but consists of noise or 
other incidental spurious signals, no remarkable correlation peaks should 
arise. The code phase difference is seen from the row index of this highest 
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value, and correspondingly the frequency deviation is seen from the column 
index. The example of Fig. 6 illustrates a situation in which a value, which is 
clearly higher than the other values, has been found. If no such value is 
found in the incoherent search matrix S^, which means that a signal 

5 transmitted by the satellite searched for was probably not received in the 
frequency range examined, the frequency range to be examined Is changed, 
and the above mentioned steps are carried out for forming an incoherent 
search matrix. With this method, the whole range of 6 kHz can be examined 
by repeating the steps described above as many times as Is required. 

10 

When required, the steps described above can be repeated for the whole 
frequency range to be examined, and the incoherent search matrixes formed 
at different times of repetition or only the possible peaks can be saved 
before searching for the highest correlation peak. Then the possibility of 
15 false interpretations can be reduced In such a situation, for instance, that the 
threshold value has been set too low and a spurious signal can cause a 
false interpretation. 

If there were devices available, which could perform the code phase and 
20 frequency deviation with signals, the strength of which is in the range of 
-150 dBm, the frequency range should be scanned In portions of tens of 
hertzes in width, before even an attempt could be made to determine the 
phase and the frequency deviation. With the present invention, it is possible 
to scan the frequency range in portions of hundreds of hertzes in width, 
25 which is much faster than with prior art receivers. 

After the right frequency deviation and code phase have been determined, 
the receiver can be set to the tracking mode. Reception of data is not 
successful with weaker signals, and data received via the mobile phone 

30 network, for example, must be used Instead, which is a situation known as 
such. Distance measurement is still possible with lower accuracy. The 
tracking mode is set by changing the switches (Fig. 1) to another position, 
whereby the received information is led to the tracking block 1 1 , in which a 
feedback is also formed for the fine adjustment of the frequency of the 

35 numerically controlled oscillator 5. 

For calculating the location, the receiver performs the signal reception 
preferably on the basis of signals received from at least four satellites. Then 
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the synchronization described above is repeated when required for the 
signal of each satellite, whereby the code of the satellite with which 
synchronization is performed is always selected as the reference sequence 
r(x). 

Most of the blocks needed in the implementation of the method can be 
implemented in the digital signal processor (not shown), for example. Either 
hardware-based solutions or software implementations of the digital signal 
processing unit can be used in the performance of FFT transformations. In 
addition, a control element, preferably a microprocessor or the like can be 
used to control the operation of the receiver. 

An electronic device 24 according to a preferred embodiment of the 
invention, which comprises a mobile station and a location determination 
receiver, is shown in Fig. 7. The first antenna 13 is used for the reception of 
the signal transmitted by the positioning satellites. The received signal is led 
to the first radio part 14, in which the signal is down-converted and digitized. 
The first radio part comprises, among other things, the transformation block 
2, digitizing block 3 and multiplier block 4 of the receiver shown in Fig. 1 . 
The digitized signal, which at this stage preferably comprises the I and Q 
components, is led to the digital signal processing unit 15, in which, among 
other things, sample vectors are formed in the sample vector formation block 
12. Samples are saved in the first memory media 16, which comprise, for 
instance, read/write memory and preferably also read memory and/or non- 
volatile read/write memory for saving the program code of the digital signal 
processing unit 15. In this embodiment, the signal processing unit 15 also 
comprises operations of the synchronization block 6, for instance, such as 
the formation of the correlation function matrix Cx,k with, for example, 

Fourier transformators FFT1, FFT2, .... FFTN and/or matched filters. The 
digital signal processing unit 15 preferably also includes the formation of a 
coherent search matrix Ax,k and the performance of an incoherent summing 

step. The digital signal processing unit 15 transmits information about the 
calculated phase difference and the frequency deviation to the processor 
block 1 7, which comprises, for example, a microprocessor and I/O logic. The 
processor block 1 7 controls the scanning block 8 and the first switch 9. The 
tracking block 11 is preferably at least partly implemented as program 
instructions of the processor block. The second memory media 1 8 are used 
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as the data memory and program memory of the processor block 17. It is 
obvious that the first memory media 16 and the second memory media 18 
can also comprise common memory. Location determination information can 
be shown to the user on the display 1 9. 

5 

Operations of the mobile station have also been implemented in the 
application software of the processor block 17. The display 19 can then be 
used for showing telephone call information in a manner known as such. 
With the keyboard 20, the user can control the location determination 
10 receiver and the mobile station. The coding and decoding of audio signals is 
performed by means of a codec 21. The radio part 22 and the second 
antenna 23 of the mobile station are also shown in Fig. 7. 

The present invention is not limited to the above described embodiments 
15 only, but its details can be modified without departing from the scope defined 
by the attached claims. 
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COaiims 

1. A method for synchronizing a receiver (1) with the transmitted code- 
modulated spread spectrum signal, in which method at least one reference 
code (r(x)) is used, which reference code corresponds to a code used in the 
modulation, the frequency shift of the transmitted signal and the code phase 
of the code used in the modulation is determined, dharacttorized in that the 
method includes at least the following steps: 

- sample vector formation step, in which samples (101) are taken from the 
received signal for forming the sample vectors (Pk(1), PK(2)...pk(N)), 

- correlation step (102. 103, 104, 105) for forming a correlation function 
matrix (Cx.k ) on the basis of the sample vectors (Pk(1), Pk(2)...pk(N)), 

- analysis step (106, 107) for forming a coherent search matrix (A^.k) on 
the basis of the transpose of the correlation function matrix (Cx,k ). and 

- determination step (110), in which the frequency shift and code phase are 
determined by using the values {a^,k{iJ)) of the elements of said coherent 
search matrix ( A^.k )• 

2. A method according to Claim 1 , ©Ihairacll©ir5zedl in that in the method the 
frequency range to be examined is specified, the specified frequency range 
is divided into two or more parts, whereby said sample vector formation step, 
correlation step and analysis step are performed on each part, values of the 
elements ia^,k{iJ)) of the coherent search matrix (A^.k) are saved, 

whereby the determination step Is performed after the examination of the 
specified frequency range, and the frequency shift and code phase are 
determined on the basis of the highest value. 

3. A method according to Claim 1 or 2, clhiairaclteirDzedl in that a threshold 
value is specified in the method, whereby quantity values of the 
elements (a^jt(i,j)) of the coherent search matrix (A^.k) that exceed said 

threshold value are used in the determination step for determining the 
frequency shift and the code phase. 

4. A method according to any one of the Claims 1 , 2 or 3, characJeirijzed 
in that in the method said sample vector formation step, correlation step and 
analysis step are repeated (108) for forming at least two coherent search 
matrixes (A^.k), and that a summing step (109) is also performed in the 
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method, in which summing step an incoherent search matrix (s^) is formed 
by summing incoherently the values (a^^ /^(/.y)) of the equivalent elements of 
the coherent search matrix ( Ax,k ) formed at each time of repetition, whereby 
said incoherent search matrix (s^) is used in said determination step for 
5 determining the frequency shift and the code phase. 

5. A method according to Claim 4, characlieirSzeca in that in the method the 
frequency range to be examined Is specified, the specified frequency range 
is divided into two or more parts, whereby said sample vector formation step, 

10 correlation step, analysis step and summing step are performed on each 
part, and values of the elements of the incoherent search matrix (s^) are 

saved, whereby the determination step is performed after the examination of 
the specified frequency range, and the frequency shift and code phase are 
determined on the basis of the highest value. 

15 

6. A method according to Claim 4 or 5, clhiairac«eir5z©dl in that a threshold 
value is determined in the method, whereby quantity values of the 
elements of the incoherent search matrix (s^) that exceed said threshold 

value are used in the determination step for determining the frequency shift 
20 and the code phase. 

7. A method according to any one of the Claims 1 to 6, clhairaclleirSzedl in 
that the correlation function matrix (Cx,k) is formed with matched filters 

(26, 27, 28). 

25 

8. A method according to any one of the Claims 1 to 7, cBiairacterizedl In 
that 

- in the correlation step, a time-to-frequency transformation (R(x)) of a 
inverse code corresponding to said reference code (r(x)), and a time-to- 

30 frequency transformation (Pk(i)) of each sample vector 
(Pk(1), Pk(2)...pk(N)) is performed, a multiplication between the time-to- 
frequency transformation (R(x)) of said inverse code and the time-to- 
frequency transformation (Pk(i)) of each sample vector 
(Pk(1), Pk(2)...pk(N)) is performed, a frequency-to-time transformation Is 

35 performed on each multiplication result (i\/lx,k(i)), and a correlation function 
matrix (Cx,k) is formed of the frequency-to-time transformations (mx.k(i)) 

of the multiplication results (Mx,k(i)), and that 
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- a coherent search matrix (Ax,k) is formed in the analysis step by 

performing a time-to-frequency transformation on the rows of the 
transpose of said correlation function matrix (Cx,k ) • 

9. A method according to any one of the Claims 1 to 7, clhiaiiractteiriized in 
that 

- in the correlation step, a time-to-frequency transformation (R(x)) of said 
reference code (r(x)), and a time-to-frequency transformation (^(/)) of an 

inverse code corresponding to each sample vector (Pk(1), Pk(2)...pk(N)) is 
performed, a multiplication between the time-to-frequency 
transformation (R(x)) of sajd reference code (r(x)) and the time-to- 
frequency transformation (f^(/)) of the inverse code of each sample 

vector (Pk(1), Pk(2)...pk(N)) is performed, a frequency-to-time 
transformation Is performed on each multiplication result (i\/ix,k(i)), and a 
correlation function matrix (Cx.k) is formed of the frequency-to-time 

transformations (mx.k(i)) of the multiplication results (i\/ix.k(i)), and that 

- a coherent search matrix (Ax.k) is formed In the analysis step by 

performing a time-to-frequency transformation on the rows of the 
transpose of said correlation function matrix (Cx,k) • 

10. A method according to any one of the Claims 1 to 7, dhiaracHerized In 
that 

- in the correlation step, a time-to-frequency transformation (R(x)) of said 
reference code (r(x)), and a time-to-frequency transformation (Pk(/)) of 
each sample vector (Pk(1), Pk(2)...pk(N)) is formed, a complex 
conjugate (R*(x)) is formed of the time-to-frequency transformation (R(x)) 
of the reference code (r(x)), a multiplication between said complex 
conjugate and the time-to-frequency transformation (Pk(/)) of each 
sample vector (Pk(1). Pk(2)...pk(N)) is performed, a frequency-to-time 
transformation is performed on each multiplication result (Mx,k(i)), and a 
correlation function matrix (Cx,k) is formed of the frequency-to-time 

transformations (mx.k(i)) of the multiplication results (Mx.k(i)), and that 

- a coherent search matrix (Ax.k) is formed in the analysis step by 

performing a time-to-frequency transformation on the rows of the 
transpose of said correlation function matrix (Cx.k) • 
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11. A method according to any one of the Claims 1 to 7, clhairacterDzed in 
that 

- in the correlation step, a time-to-frequency transformation (R(x)) of said 
reference code (r(x)), and a time-to-frequency transformation (Pk(/)) of 
each sample vector (Pk(1), Pk(2)...pk(N)) is formed, a complex 
conjugate (P*k(i)) is formed of the time-to-frequency 
transformation (Pk(/))) of each sample vector (pk(1), Pk(2)...pk(N)), a 
multiplication between each said complex conjugate and time-to- 
frequency transformation (R(x)) of the reference code (r(x)), is performed, 
a frequency-to-time transformation Is performed on each multiplication 
result (Mx,k(i)), and a correlation function matrix (Cx,k) is formed of the 

frequency-to-time transformations (mx.k(i)) of the multiplication results 
(i^x.k(i)). and that 

- a coherent search matrix (Ax.k) is formed in the analysis step by 

performing a time-to-frequency transformation on the rows of the 
transpose of said correlation function matrix (Cx.k) - 

12. A method according to any one of the Claims 8 to 1 1 , clhairacterozed in 
that said time-to-frequency transformation is a Fourier transformation, 
whereby said frequency-to-time transformation is an inverse Fourier 
transformation. 

13. A method according to any one of the claims 1 to 12, characHeiriized in 
that the quantity values of the elements (axk{i,j)) of the coherent search 
matrix (Ax,k) are used In said determination step (110). 

14. A method according to Claim 13, clhiairacllerDzed in that the quantity 
values of the elements (a^^ki^J)) of the coherent search matrix (Ax.k) 

raised to the second power are used in said determination step (110). 

15. A method according to any one of the Claims 1 to 14, in which the 
location of the receiver (1) is determined and the location information is 
saved, clhiairaclteriized in that the frequency range to be examined is 
specified in the method, and a starting frequency for starting the examination 
of the frequency shift is selected from said frequency range by using the 
location information saved in the receiver (1). 
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16. A method according to any one of the Claims 1 to 15, in which method 
the location and location information of the receiver (1) is determined from a 
spread spectrum modulated signal transmitted by two or more satellites, and 
In which method a unique code Is selected for each of said satellites, 

5 clhairactteirozedl in that the method includes carrying out said sample vector 
formation step, correlation step, analysis step and determination step on the 
signal transmitted by each of said satellites. 

17. A receiver (1), which comprises at least synchronization means for 
10 synchronizing the receiver (1) with the transmitted code-modulated spread 

spectrum signal, which receiver (1) includes means (16) for using at least 
one reference code (r(x)) In connection with the synchronization, which 
reference code (r(x)) corresponds to a code used in the modulation, and 
means (15) for determining the frequency shift and the code phase of the 
15 code used in the modulation, characterDzed in that the receiver (1) also 
comprises: 

- sample vector formation means (12, 16) for forming sample vectors 
(Pk(1), Pk(2)...pk(N)) from the received signal, 

- correlation means (6, FFT1, FFT2...FFTN) for forming a correlation 
20 function matrix (Cx.k) on the basis of the sample vectors 

(Pk(1), Pk(2)...pk(N)), 

- analysis means (6, 12) for forming a coherent search matrix ( Ax.k ) on the 
basis of the correlation function matrix (Cx,k ). and 

- determination means (7, 17) for determining the frequency shift and code 
25 phase by using the values of the elements (a;^ ;t of said coherent 

search matrix (Ax,k)- 

18. A receiver (1) according to Claim 17, cliiairac«©irii2©d in that it comprises 
means (17) for specifying the frequency range to be examined, means (5, 8, 

30 9) for dividing the specified frequency range Into two or more parts, whereby 
the formation of sample vectors, the formation of a correlation function 
matrix (Cx.k) and the formation of a coherent search matrix (Ax,k) are 

arranged to be performed on each part, means (16) for saving the values of 
the elements (a;, ;t('.y)) of the coherent search matrix (Ax,k). and means 

35 (15) for determining the frequency shift and code phase on the basis of the 
highest value. 
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19. A receiver (1) according to Claim 17 or 18,. clhiairaclteiriizedl in that it 
comprises means (15) for specifying the threshold value, and means (15) for 
comparing the values of said threshold value and the values of said 
coherent search matrix (Ax.k) for determining the frequency shift and code 

5 phase. 

20. A receiver (1) according to any one of the Claims 17, 18 or 19, 
chairaclleirDzed in that the formation of sample vectors, the formation of a 
correlation function matrix (Cx,k) and the formation of a coherent search 

10 matrix (Ax,k) are arranged to be repeated at least two times for forming a 
coherent search matrix (A^.k), and that the receiver (1) also comprises 
summing means for forming an incoherent search matrix (s^) by summing 
the values (a;f /f(/,y)) of the equivalent elements of the coherent search 
matrix (A^.k) formed at each time of repetition, whereby said Incoherent 

15 search matrix (5^^) is used in said determination step for determining the 
frequency shift and code phase. 

21 . A receiver (1) according to Claim 17, clhairaclteirozed in that it comprises 
means (17) for specifying the frequency range to be examined, means (5, 8, 

20 9) for dividing the specified frequency range into two or more parts, whereby 
the formation of sample vectors, the formation of a correlation function 
matrix (Cx,k) and the formation of a coherent search matrix (Ax,k) are 

arranged to be performed on each part, and summing of the coherent 
search matrix (Ax,k) to the incoherent search matrix (s^), means (16) for 

25 saving the values of the elements of the incoherent search matrix {s^), and 
means (15) for determining the frequency shift and code phase on the basis 
of the highest value. 

22. A receiver (1) according to Claim 21 , clhiairacfteirozed in that it comprises 
30 means (1 5) for specifying the threshold value, and means (1 5) for comparing 

the values of said threshold value and the values of said incoherent search 
matrix ( ) for determining the frequency shift and code phase. 

23. A receiver (1) according to any one of the Claims 17 to 22, 
35 characteiTDzed in that it comprises matched filters (26, 27, 28) for forming a 

correlation function matrix (Cx,k )■ 
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24. A receiver (1) according to any one of the Claims 17 to 23, 
clhiairacteiriized in that 

- the correlation means cojTiprise means (16) for forming a time-to- 
frequency transformation ( R(x)) of an inverse code corresponding to said 

reference code (r(x)), and a time-to-frequency transformation (Pk(i)) of 
each sample vector (Pk(1), Pk(2)...pk(N)), means (15) for performing a 
multiplication between the time-to-frequency transformation ( R(x)) of said 

inverse code and the time-to-frequency transformation (Pk(i)) of each 
sample vector (Pk(1), Pk(2)...pk(N)), means (15) for performing a 
frequency-to-time transformation on each multiplication result (i\/ix,k(i)), 
and means (15) for forming a correlation function matrix (Cx.k) of the 

frequency-to-time transformations (mx,k(i)) of the multiplication results 
(Mx,k(i)), and that 

- the analysis means comprise means (15) for forming a coherent search 
matrix(Ax.k) by performing a time-to-frequency transformation on the 
rows of the transpose of said correlation function matrix (Cx,k )- 

25. A receiver (1) according to any one of the Claims 17 to 23, 

clhairactteirDzecil in that 

- the correlation means comprise means (16) for forming a time-to- 
frequency transformation (R(x)) of said reference code (r(x)) and a time- 
to-frequency transformation iFU{i)) of an Inverse code corresponding to 

each sample vector (Pk(1), Pk(2)...pk(N)), means (15) for performing a 
multiplication between the time-to-frequency transformation (R(x)) of said 
reference code (r(x)) and the time-to-frequency transformation (^(/)) of 

the inverse code of each sample vector (Pk(1), Pk(2)...pk(N)), means (15) 
for performing a frequency-to-tlme transformation on each multiplication 
result (i\/ix,k(i)), and means (15) for forming a correlation function 
matrix (Cx,k) of the frequency-to-time transformations (mx.k(i)) of the 

multiplication results (i\/ix.k(i)), and that 

- the analysis means comprise means (15) for forming a coherent search 
matrix(Ax,k) by performing a time-to-frequency transformation on the 
rows of the transpose of said correlation function matrix (C^.k )- 

26. A receiver (1) according to any one of the Claims 17 to 23, 
characterized in that 
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- the correlation means comprise means (16) for forming a complex 
conjugate (R*(x)) of the time-to-frequency transformation (R(x)) of the 
reference code (r(x)) and a time-to-frequency transformation (^(/)) of 

each sample vector (Pk(1), Pk(2)...pk(N)), means (15) for performing a 
5 multiplication between the complex conjugate (R*(x)) of the time-to- 
frequency transformation of the reference code (r(x)) and the time-to- 
frequency transformation (^(/)) of each sample vector 

(Pk(1). Pk(2)...pk(N)), means (15) for performing a frequency-to-tlme 
transformation on each multiplication result (Mx.k(i)), and means (15) for 
10 forming a correlation function matrix (Cx,k) oi the frequency-to-time 

transformations (mx,k(i)) of the multiplication results (l\/lx.k(i)), and that 

- the analysis means comprise means (15) for forming a coherent search 
matrix(Ax,k) by performing a time-to-frequency transformation on the 
rows of the transpose of said correlation function matrix (Cx.k )■ 

15 

27. A receiver (1) according to any one of the Claims 17 to 23, 

chairaclteD'azGd in that 

- the correlation means comprise means (16) for forming a complex 
conjugate (P*k(i)) of the time-to-frequency transformation (^(/))of each 

20 sample vector (Pk(1 ), Pk(2)...pk(N)) and a time-to-frequency transformation 
(R(x)) of the reference code (r(x)), means (15) for performing a 
multiplication between the time-to-frequency transformation (R(x)) of the 
reference code (r(x)) and the complex conjugate (P*k(i)) of the time-to- 
frequency transformation of each sample vector (Pk(1), Pk(2)...pk(N)), 

25 means (15) for performing a frequency-to-time transformation on each 
multiplication result (Mx,k(i)), and means (15) for forming a correlation 
function matrix (Cx,k) of the frequency-to-time transformations (mx.k(i)) of 

the multiplication results (Mx.k(i)), and that 

- the analysis means comprise means (15) for forming a coherent search 
30 matrix(Ax.k) by performing a time-to-frequency transformation on the 

rows of the transpose of said correlation function matrix (Cx.k )■ 

28. A receiver (1) according to any one of the Claims 24 to 27, 
clhairaclteirizedl in that said time-to-frequency transformation is a Fourier 

35 transformation, whereby said frequency-to-time transformation is an inverse 
Fourier transformation. 
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29. A receiver (1) according to any one of the Claims 17 to 28, 
characllerilzecJ in tliat quantity values, such as absolute values, of the 
elements (flj,^^ (/,;•)) of the coherent search matrix (Ax,k) are used in the 

determination means. 

30. A receiver (1) according to any one of the Claims 17 to 29, 
©haradtenzed in that quantity values of the elements (a^jt(i,y)) of the 
coherent search matrix (Ax,k) raised to the second power are used in the 
determination means. 

31. A receiver (1) according to any one of the Claims 17 to 30, which 
comprises means (14, 15, 17) for determining the location of the receiver (1) 
and means (18) for saving the location information, cUiairacireirozeidl in that 
the receiver (1) also comprises means (17) for specifying the frequency 
range to be examined, and means (17, 22, 23) for selecting the starting 
frequency from said frequency range on the basis of the location information 
saved in the receiver (1). 

32. An electronic device (24), which comprises at least a location 
determination receiver (1), synchronization means for synchronizing the 
location determination receiver (1) to the transmitted code-modulated spread 
spectrum signal, which location determination receiver (1) includes means 
(16) for using at least one reference code (r(x)) in connection with the 
synchronization, which reference code (r(x)) corresponds to a code used in 
the modulation, and which electronic device comprises means (15) for 
determining the frequency shift of the transmitted signal and the code phase 
of the code used in the modulation, chairaclleirizedl in that the electronic 
device (24) also comprises: 

- sample vector formation means (12, 16) for forming sample vectors 
(Pk(1), Pk(2)...pK(N)) from the received signal, 

-correlation means (6, FFT1, FFT2...FFTN) for forming a correlation 
function matrix (Cx,k) on the basis of the sample vectors 

(Pk(1), Pk(2)...pk(N)), 

- analysis means (6, 12) for forming a coherent search matrix (Ax,k) on the 
basis of the correlation function matrix (Cx,k ). and 
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- determination means (7, 17) for determining the frequency shift and code 
phase by using said values of the elements {a^^k{i^J)) of the coherent 
search matrix ( Ax,k )■} 

5 33. An electronic device (24) according to Claim 32, comprising means (14, 
15, 17) for determining the location of the electronic device (24) and means 
(18) for saving the location information, clhaiirsiclteirDzedl in that the electronic 
device (24) also comprises means (17) for specifying the frequency range to 
be examined, and means (17, 22, 23) for selecting the starting frequency 
10 from said frequency range on the basis of the location information saved in 
the receiver (1). 

34. An electronic device (24) according to the Claim 32 or 33, 
cDnairsictteirDzedl in that it also comprises means for performing mobile station 

15 operations. 

35. An electronic device (24) according to Claim 34, clhairacierozedl in that it 
comprises means (22, 23) for establishing a data transfer connection to a 
mobile communication network, whereby the means for determining the 

20 location of the electronic device (24) comprise means (22. 23) for retrieving 
information used in the location determination from the mobile 
communication network. 
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(57) Abstiractt: 

The invention relates to a metliod for synchronizing a 
receiver (1) with a transmitted code-modulated spread 
spectrum signal. At least one reference code (r(x)) is 
used in the method, which reference code 
corresponds to a code used in the modulation. The 
frequency shift of the transmitted signal and the code 
phase of the code used in the modulation are 
determined in the method. The method also includes 
at least the following steps: 

- sample vector formation step, in which samples 
(101) are taken from the received signal for forming 
the sample vectors (pk(1), Pk(2)...pk(N)), 

- correlation step (102, 103, 104, 105) for forming a 
correlation function matrix (Cx,k) on the basis of the 

sample vectors (Pk(1). Pk(2)...pk(N)). 

- analysis step (106, 107) for forming a coherent 
search matrix (A^.k) on the basis of the transpose of 
the correlation function matrix (Cx,k ). and 

- determination step (110), in which the frequency shift 
and code phase are determined by using the values of 
the elements (a^^kiiJ)) of said coherent search 

matrix(Ax,k)- 
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